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1 .  PLASMA  IGNITION 


The  ignition  of  a  plasma  over  a  target  has  been  found  to  consist  of 
two  steps,  the  vaporizing  of  a  small  quantity  of  the  target  material  followed 
by  the  breakdown  of  the  "vapor.  The  intensities  are  of  course  below  the  break¬ 
down  thresholds  of  clean  air.1  The  production  of  aluminum  vapor  can  arise 
irom  the  heatina  of  the  bulk  aluminum  material  and  this  is  indicated  in  the 
theory  plot  of  Fig.  1.1  by  the  line  labeled  bulk.  This  heating  process  in 
which  the  surface  temperature  is  proportional  to  the  intensity  and  the  square 
root  of  the  time  T  *  i/t  when  combined  with  the  fact  that  $  =  It  for  top  hat 
pulses,  leads  to  a  criterion  that  $1  >  6x10®j2/cm^  s.  The  other  mechanism  for 
producing  vapor  that  has  been  identified  for  aluminum  is  the  vaporization  of 
flakes,  small  thermally  insulated  pieces  of  aluminum  lying  on  the  aluminum 
surface.  It  has  been  found  that  2  J/cm2  is  required  to  vaporize  such  small 
flakes,1  (independent  of  intensity)  and  this  point  appears  as  the  horizontal 
line  labeled  FLAKES  in  Fig.  1.1.  The  curves  marked  1 -D  and  3-D  refer  to  the 
breakdown  in  the  vapor  where  it  is  seen  that  the  time  to  breakdown  the  vapor 
is  small  compared  to  the  time  to  generate  the  vapor  at  high  intensities  and 
can  be  neglected. 

When  the  breakdown  threshold  curves  are  compared  to  the  conditions 
obtainable  in  the  present  10.6  y  tests.  Fig.  1.2,  it  is  seen  that  the  flake 
limit  should  prevail.  There  was  considerable  scatter  in  the  results  as  an 
exact  determination  of  the  ignition  time  was  quite  difficult,  however  the  data 
did  tend  to  center  around  a  value  of  2  J/cm2,  Fig.  1.3,  indicating  that 
flakes  did  indeed  dominate  the  ignition  process. 

The  situation  for  grafoil  is  less  clear.  The  ignition  times  appear 
very  similar  to  those  on  aluminum  with  an  ignition  on  fluence  of  perhaps 
4  J/cm2  being  more  appropriate.  Using  a  philosophy  similar  to  that  for  alumi¬ 
num,  namely  that  the  ignition  process  is  dominated  by  the  production  of  the 
vapor,  this  value  is  quite  consistent  with  our  estimates  of  the  joules  per  cm2 
required  to  produce  carbon  vapor.  However  predictions  of  the  time,  and  there¬ 
fore  energy,  required  to  breakdown  the  vapor  were  not  negligible  as  in  the 
aluminum  case,  and  the  physical  mechanisms  responsible  for  the  breakdown  of 
the  carbon  vapor  must  be  considered  unknown  at  the  present  time.  The  details 
of  PSI's  examination  of  the  problem  are  covered  in  the  next  section. 
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1.1 


Laser  Induced  Breakdown  of  Grafoil"  Surfaces 


In  this  section,  we  predict  the  threshold  intensity  (or  fluence)  to 
cause  breakdown  above  a  Grafoil”  surface.  The  wavelength  considered  is  10.6 
ym  and  the  pulse  shape,  corresponding  to  the  CO2  laser  used  at  NRL,  is  shown 
in  Fig.  1.4. 


Our  model  for  breakdown  is  as  follows.  The  laser  energy  absorbed  in  a 
small  layer  of  Grafoil”  material  below  the  surface  results  in  a  heating  of 
this  layer  till  vaporization  occurs.  The  vapor  that  comes  off  the  surface 
drives  a  shock  into  the  air.  We  assume  that  breakdown  occurs  in  the  carbon 
vapor  due  to  inverse  bremsstrahlung  absorption  by  the  electrons  in  the  vapor. 
This  causes  an  increase  in  vapor  temperature,  which  in  turn  causes  an  increase 
in  the  number  density  of  electrons.  Breakdown  is  assumed  to  have  occurred  if, 
by  the  end  of  the  pulse,  a  critical  ionization  level  has  been  achieved  (ne/n  > 
1%).  The  work  performed  here  is  an  extension  to  10.6  ym  of  a  model  developed 
for  shorter  wavelengths  (X  =  1.06  ym  and  X  =  0.35  ym).2  In  Section  1.1,  we 
analyze  the  heating  of  the  material  to  vaporization.  We  then  discuss  in  Sec¬ 
tion  1.2,  the  vapor  properties  before  breakdown.  Our  vapor  breakdown  model  is 
presented  in  Section  1.3  and  results  of  specific  calculations  for  the  pulse 
shape  shown  in  Fig.  1.4  are  presented  in  Section  1.4. 


1 .2  Heating  to  Vaporization 

The  response  of  the  solid  material  to  laser  radiation,  in  the  absence 
of  radiation  and  vaporization  losses  is  given  by  the  heat  conduction  equation 


C  3T 
3t 


ill 

3xJ 


(1  -  R)  -x/l 


Ie 


(1.1) 


where  p,  C  and  K  are  the  density,  specific  heat  and  thermal  conductivity  of 
the  material,  I  is  the  incident  laser  intensity,  T  is  the  temperature,  R  is 
the  surface  reflectivity  and  l  the  absorption  length  in  the  material. 

Grafoil”  is  a  form  of  graphite  in  which  thin  lamellae  of  graphite 
layers  composed  of  microscopic  crystals  all  having  their  c  axis  perpendicular 
to  the  lamellae  are  pressure  bonded  to  each  other.  Though  pure  graphite  has  a 
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density  of  2.2  g/cm3,  the  density  of  Grafoil”  is  a  factor  of  two  lower  due  to 


the  voids  between  lamellae.  The  thermal  properties  of  Grafoil”  are  shown  in 
Table  1.1.  We  have  also  shown  in  this  table  the  thermal  properties  of  a  pure 
graphite  single  crystal,  as  well  as  the  properties  of  commercial  grade  poly¬ 
crystalline  ATJ  Graphite,  which  is  actually  a  composite  of  micro  crystals  with 
preferred  orientation  due  to  the  forming  process.  Note  that  a  single  graphite 


crystal  has  a  thermal  conductivity  in  the  a  direction  that  is  one  hundred 


times  that  in  the  c  direction  and  that  is  comparable  to  the  thermal  conduc¬ 
tivity  of  a  metal.  Thermal  conductivity  of  bulk  graphite,  which  is  composed 
of  many  crystals,  is  less  anisotropic.  The  conductivity  of  Grafoil”  in  the 
direction  normal  to  the  lamellae,  is  much  smaller  than  that  of  monocrystalline 


graphite  in  the  c  direction  due  to  the  voids  between  lamellae. 


A  simplification  of  Eq.  (1.1)  can  be  obtained  by  comparing  the  thermal 
diffusion  length  ZD  =  [k  Tp/tpc)]1/^  during  the  pulse  time  (tp)  to  the  absorp¬ 
tion  length  l.  If  ZD  <<  l  then  thermal  conduction  during  the  pulse  is  negli¬ 
gible  and  the  temperature  at  any  time  during  the  pulse,  obtained  by  integra¬ 


tion  of  Eq.  (1.1),  is: 


/  p ( T)C( T)dT  = 


(1  -  R)  -x/Z 
— 1 -  6 


/  I  dt' 


h(T( x) )  =  h(To)  +  ^  R  e~X/£  *(t) 


( 1  .2) 


where  h  is  the  enthalpy  per  unit  mass,  $  the  incident  fluence  and  when  we 
neglected  variation  of  absorption  length  and  reflectivity  with  temperature. 
Equation  (1.2)  can  readily  be  solved  to  obtain  T  as  a  function  of  4>.  We  show 


in  Fig.  1.5  the  complex  index  of  refraction  of  graphite  as  measured  by  Foster 
and  Howarth.3  Extrapolating  the  curves  up  to  X  =  10.6  ym,  we  find  n  =  n-|+in2 


=  3.9+i  0.7  from  which  we  calculate  the  reflection  coefficient 


R  -  -U-l ."i  ,  0.36 
(n,  +  1)  +  n  2 


TABLE  1 . 1 

Thermal  Properties  of  Graphite 


Material 

P 

( g/cm3 ) 

C 

(J/gK) 

Kc 

W/(cm  K) 

Ka 

W/(cm  K) 

Grafoil+ 

1.1 

1  .8 

0.03 

0.42 

Polycrystalline 

ATJt  Graphite 

2.2 

1  .8 

0.3 

0.4 

Graphite  Single* 
Crystal 

2.2 

0.7  (T=300  K) 

2  (T  =  1000°C) 

0.08 

20.0 

"•Union  Carbide  Data 

Sheet, 

T  =  1100°C 

tTouloukian,  page  1,  T  =  1100°C 

*A.  J.  deCombarieu,  J.  Phys.  (Prance)  28,  931  (1968) 


3 


-  I 

tfl 


e 


5 


i  ! 


£ 


'*  £ 


I 


5  | 
.-v  3 

4  5 


-H 


and  the  absorption  length 


l  =  X/(4nn2)  *  1.2  in 

This  length  is  to  be  compared  with  the  diffusion  distance  during  the  pulse 
Taking  Tp  =  2  *  10“^  s  and  considering  a  monocrystal  near  the  surface,  we  find 
in  the  c  direction  iD  =  0.14  urn  and  in  the  a  direction  £D  =  1.2  pm.  If  the 
lateral  (i.e.,  in  the  a  direction  parallel  to  the  surface)  dimensions  of  the 
micro  crystals  exceed  a  few  microns  the  response  of  Grafoil™  during  the  laser 
pulse  will  be  the  same  as  that  of  an  aggregate  of  microcrystals  that  are  un¬ 
coupled  from  each  other.  The  main  point  to  be  made  here,  however,  is  that 

whatever  thermal  conductivity  we  choose  in  the  c  direction  (normal  to  the  sur¬ 

face),  we  will  have  iD  <<  l  so  that  Eq.  (1.2)  holds.  The  surface  temperature 
reached  by  the  Grafoil™  is  shown  in  Fig.  1.6  as  a  function  of  incident  fluence, 
using  i  =  1.2  ym,  R  =  0.4  and  tabulated  values  of  C  for  graphite.  An  incident 

fluence  of  3.3  J/cm^  is  required  to  heat  the  surface  to  T  =  4000  K  at  which 

point  the  vapor  pressure  reaches  1  atm. 

1 .3  Vaporization  of  Grafoil1* 

Once  the  vapor  pressure  above  the  surface  exceeds  1  atm  (i.e.,  Ts  > 

4000  K) ,  the  carbon  vapor  that  leaves  the  surface  drives  a  shock  into  the  air. 
We  will  have  a  situation  as  shown  in  Fig.  1.7.  Zone  1  is  the  region  of  shocked 
air.  Zone  2  is  the  region  where  the  vapor  is  expanding  at  a  uniform  velocity 
u  away  from  the  surface  and  is  in  pressure  equilibrium  with  the  shocked  air. 
Zone  3  is  the  Knudsen  layer  where  the  half  Maxwellian  distribution  in  velocity 
space  of  particles  leaving  the  surface  is  converted  into  a  full  Maxwellian  cen¬ 
tered  around  u.  We  use  the  model  developed  by  C.  Knight^'6  to  solve  for  the 
mass  flow  rate  as  a  function  of  surface  temperature  Ts.  The  initial  composi¬ 
tion  of  the  vapor  as  a  functon  of  Ts  is  assumed  to  be  the  same  as  that  of  sat¬ 
urated  vapor  at  Tg  and  is  obtained  by  use  of  a  chemical  equilibrium  code  that 
uses  the  JANAF  thermochemical  properties  of  the  vapor  constituents.  We  show 
the  concentration  of  the  vapor  species  in  Fig.  1.8.  The  electron  concentration 
is  given  by  the  balance  of  the  two  reactions: 

C  t  C+  +  e  e  *  11.3  eV  (1 .3a) 

e  +  C  t  C  ~  t  =  4  eV  (1  .3b) 

&  c  3 
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LAYER  SURFACE 


Fig.  1.7  Geometry  for  vapor  expansion. 


n®,  the  electron  density,  is  well  approximated  over  the  temperature  range 
studied  by  the  relation1 


ne  =  3.2  x  1016  T3/2  exp( -92 ,800/T) 


(1  .4) 


Hie  vapor  blowing  rate  m  was  obtained  by  solving  the  Knudsen  layer 
and  continuity  equations  given  by  C.  Knight. ^  a  computer  program  that  solves 
ttiese  equations  was  developed  at  PSI  by  Les  Popper.2  We  show  the  results  in 
Pig.  1.9.  The  y  of  the  vapor  was  taken  to  be  1.2  and  the  saturated  vapor 
pressure  was  taken  to  be 

p(Tg)  =  9  x  io9  exp-(94,000/Tg)  atm  . 

The  vapor  chokes  and  the  flow  becomes  supersonic  when  the  temperature  exceeds 
5100  K. 

In  the  absence  of  laser  induced  breakdown  in  the  vapor  (causing  ab¬ 
sorption  of  the  laser  beam  by  the  vapor),  the  surface  temperature  eventually 
reaches  an  asymptotic  value  for  which  the  rate  of  enthalpy  loss  by  vaporiza¬ 
tion  balances  the  power  absorbed  from  the  laser  beam  at  the  surface.  We  have 
plotted  this  rate  of  enthalpy  loss  by  multiplying  the  mass  loss  rate  m,  shown 
in  Fig.  1.9,  by  (A  hv  +  /  cdT' )  where  A  hv  is  the  heat  of  vaporization  (~23 
kJ/g)  and  / CdT'  the  enthalpy  rise  in  the  solid  to  reach  the  vaporization  tem¬ 
perature.  This  power  loss  rate  is  shown  in  Fig.  1.10.  We  have  also  plotted 
in  Fig.  1.10  the  radiative  losses  from  the  surface  (which  are  several  orders 
of  magnitude  lower) . 


Breakdown  of  the  Vapor 


The  electrons  absorb  laser  radiation  due  to  inverse  bremsstrahlung 
collisions  with  the  neutrals.  At  sufficiently  high  laser  fluence  some  elec¬ 
trons  may  gain  a  sufficient  amount  of  energy  that  upon  impact  they  can  ionize 
the  vapor  leading  to  avalanche  ionization.  A  proper  calculation  of  the  ioni¬ 
zation  rate  would  involve  numerical  solution  of  the  Boltzman  equation.  The 
cross  sections  for  the  various  inelastic  processes  in  carbon  are  however  un¬ 
known.  Because  the  saturated  vapor  above  a  carbon  target  is  mainly  composed 
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Fig.  1.9  Mass  flow  rate  versus  surface  temperature. 
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Fig.  1.10  Loss  rate  by  vaporization  and  radiation 
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of  molecules  and  because  the  effective  ionization  potential  of  carbon  (tj  + 
ea  =  15.3  eV)  is  close  to  that  of  ^(Ej  =  15.9  eV)  we  believe  that  laser  in¬ 
duced  breakdown  of  carbon  will  not  differ  too  significantly  from  that  of  ni¬ 
trogen  and  oxygen,  the  breakdown  thresholds  of  which  are  known.  Even  though 
the  vapor  is  hot,  super-elastic  collisions  will  not  lower  the  threshold  signi¬ 
ficantly  below  room  temperature  values  since  kT  <<  £j,  ea.  The  effect  of 
super-elastic  collisions  in  N2  has  been  studied  by  making  electron  distribu¬ 
tion  predicts  5000  K,  taking  into  account  super-elastic  collisions.  We  found 
the  threshold  at  T  =  5,000  K  to  differ  insignificantly  from  the  room  tempera¬ 
ture  value  of  3  x  109  W/cm2 . 

Based  on  the  above  arguments,  we  model  the  breakdown  of  the  vapor  at 
laser  fluxes  below  109  W/cm2  as  a  heating  of  the  vapor  through  inelastic  and 
elastic  collisions  of  vapor  molacules/atoras  with  electrons.  The  electrons,  as 
we  mentioned  above,  get  heated  by  inverse  bremsstrahlung  absorption.  The  ab¬ 
sorption  coefficient  for  collisions  with  carbon  atoms  has  been  calculated  by 
Geltman7  and  we  show  his  results  in  Fig.  1.8.  Due  to  the  large  number  of 
lower  lying  electronic  states  of  carbon  and  to  energy  losses  in  excitation  of 
rotational  and  vibrational  degrees  of  freedom  of  the  gas,  we  consider  that  the 
electron  temperature  does  not  depart  significantly  from  the  vapor  temperature. 
This  assumption  should  be  well  verified  at  fluxes  below  108  W/cm2.  We  con¬ 
sider  that  the  electron  concentration  is  initially  that  given  by  using  Saha 
equilibrium  in  the  saturated  vapor.  We  can  neglect  recombination  of  electrons 
with  ions  during  the  expansion  through  the  Knudsen  layer  when  the  electron 
density  is  below  lO1^  cm~3  (i.e.,  Ts  <  5000  K)  because  the  three  body  recombi¬ 
nation  rate  is  small.  At  higher  densities,  we  will  also  neglect  recombination 
because  the  laser  beam  will  always  maintain  Te  >  Tg.  The  heating  of  the  gas 
is  given  by  the  following  relation. 

nmcv^-£=x(1  -  e  hvATj  ng  n  j(  1  +  r)  (1.5) 

where  m  is  the  average  molecular  weight.  The  absorption  cross  section  (cor- 

— h  /kT 

rected  for  stimulated  emission)  x(1  -  e  V  )  is  plotted  in  Fig.  1.11.  We 
assume  that  the  electron  concentration  is  given  by  the  Saha  equation 


!!*I  =  AT372  e“(eI  +  ea>AT. 
n 


(1  .6) 


where  we  have  added  to  the  ionization  potential  ex  the  attachment  energy  ea 
because  of  the  second  reaction  (1.3b).  It  can  be  shown6  that  the  breakdown 
time  is  equal  to  the  time  it  takes  to  e-fold  the  electron  density  once,  assum¬ 
ing  equilibrium  gas  properties.  Taking  the  leading  term  of  the  logarithmic 
derivative  of  ne  from  Eq.  (1.6),  we  obtain 


£i  .+  .f  a 

2  kT^ 


dT 


Then  setting  d  ne/ne  =  1 ,  we  obtain  the  temperature  rise  necessary  for  one  e- 
fold  in  electron  density  and  inserting  this  into  Eq.  (1.5)  we  find 


2(kT)2 


VB 


(Y  -  1)(Cj  +  ea)X ( 1  -  e  h^/kT)  tied  +  R)I 
We  evalute  Bg.  (1.7)  using  the  following  values  of  the  parameters: 


X(1  -  e'hv/kT) 


10~37( 


5000' 


0.8 


cm5 


R  =  0.4 
Y  =  1.16 


+  ea  =  15.3 
ex  =  13  eV  = 

ne  =  n  neo  =  3.2  *  1016  rj  Ts3/2  exp( -92,800/Ts ) 


(1.7) 


where  n  is  the  density  decrease  (n  =  ne/neo)  through  the  Knudsen  layer  which 
has  the  lower  bound  of  0.3  when  the  flow  is  choked  (Ts  >  5000  K) .  We  find 


1 .4  x  10-5  f 
<J>VB  “  ItVB  =  „  l 


5000 


-0.3  r  /-50001-. 

expi  1  8.5(— —  ■  -J  }  J/cm2 
Ts 


(1  .8) 


The  threshold  flux  for  breaking  down  the  vapor  as  a  function  of  surface 
temperature  is  plotted  in  Fig.  1.12. 
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Breakdown  Prediction  for  Experiments  at  NRL 


We  are  now  in  a  position  to  make  breakdown  predictions  for  the  pulse 

i  shape  shown  in  Fig.  1.4.  The  thresholds  for  breakdown  that  we  present  in  this 

1  section  were  calculated  as  follows.  A  fluence  $i,  given  by  Eq.  (1.2),  is  re- 

(  quired  to  raise  the  surface  of  the  Grafoil"  to  Ts  and  an  extra  fluence  4>vb^ts^ 

i  given  by  Eq.  (1.8)  is  required  to  breakdown  the  vapor.  The  total  fluence  to 
breakdown  is  thus 


$ 

TOT 


+  *,ro<TS> 


VB 


(1.9) 


Given  the  pulse  profile  and  the  peak  intensity  in  the  pulse  Ip,  it  is  straight¬ 
forward  to  calculate  when  in  the  pulse  a  total  fluence  has  been  delivered 

and  therefore  to  calculate  the  time  to  breakdown.  The  fluence  <J>  in  Eq.  (1.9), 

1  however,  must  be  treated  as  a  parameter  because  during  the  pulse,  as  <t>  in¬ 

creases,  the  vapor  is  coming  off  at  higher  and  higher  pressures  and  temper¬ 
atures,  so  that  as  Ts  increases,  4>vb(ts)  decreases  (see  Fig.  1.9).  We  have 
therefore  solved  for  the  minimum  value  of  ^tot  as  $i  is  varied  in  Eq.  (1.9). 
This  gives  us  a  minimum  breakdown  time  as  a  function  of  Ip  (or  of  total  pulse 
fluence  which  is  proportional  to  Ip) .  As  Ip  is  decreased  from  some  large  val¬ 
ue,  we  find  that  breakdown  occurs  later  and  later  in  the  pulse  till  below  some 

threshold  Ith  we  cannot  find  a  solution.  Our  results  are  shown  in  Fig.  1.13, 
as  a  function  of  peak  intensity,  Ip,  in  the  pulse.  The  main  assumption  made  in 
specifying  Fig.  1.13  is  that  the  rise  in  intensity  to  I  =  Ip  occurs  fast  enough 
that  vaporization  losses  can  be  neglected  during  the  rise  time.  The  steady 
state  temperature  that  the  surface  would  reach  if  I  remained  constant  during 
the  pulse  is  shown  in  Fig.  1.13a.  The  incident  fluence,  <j)H,  required  to  heat 
the  surface  to  Ta  is  shown  in  Fig.  1.13b  and  is  in  the  4  to  6  J/cm2  over  the 
intensities  of  interest.  Hie  additional  fluence  <f>vBr  required  to  breakdown  the 
vapor  leaving  the  surface  at  Tsg  is  also  plotted  in  Fig.  1.13b.  The  sum 
<|>H  +  ‘t'VBr  becomes  less  than  the  pulse  fluence  when  Ip  >  4x1 08  W/cm2  and  this 
occurs  at  a  pulse  fluence  of  65  J/cm2. 
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Fiq.  1.13  Steady  state  surface  temperature  and  breakdown  fluence  as  a 
function  of  laser  intensity  at  focus.  (A)  Steady  state 
temperature  of  surface.  B)  fluence  to  breakdown. 


Our  theoretical  estimate  is  to  be  compared  with  an  experimentally  ob¬ 
served  threshold  of  2x10®  W/cm2  and  pulse  fluence  of  5  J/cm^.  a  factor  of  two 
difference  in  peak  intensity  may  arise  from  the  poor  measurement  of  laser  pulse 
profile  since  a  gain  switch  spike  a  factor  of  two  higher  than  the  measured  peak 
may  well  be  apparent  with  a  faster  detector.  The  larger  factor  observed  in 
impulse  fluence  to  breakdown  cannot  be  explained  in  this  manner.  We  have  used 
a  theoretical  inverse  Bremsstrahlung  absorption  cross  section  calculated  by 
Geltman  for  mona tonic  carbon.  The  main  species  in  the  vapor  is  C3  and  it  is 
not  inconceivable  that  the  cross  section  could  be  a  factor  of  3  higher  or  more. 
If  we  increase  the  cross  section  by  a  factor  of  three  we  now  calculate  for 
Ip  =  4x10®  W/cm2,  a  breakdown  fluence  <j>vBr  =  1  7  J/cm2  and  <(>H  +  <f>VBr  =  22  J/cm2 
in  much  better  agreement  with  the  data.  Note  also  that  we  gave  extended  satu¬ 
rated  vapor  pressure  data  to  temperatures  and  pressures  well  beyond  the  range 
where  the  thermodynamic  data  has  been  measured. 
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2.  PRESSURE  AT  THE  WALL 


Lindberg1  showed  that  the  damage  to  re-entry  vehicles  depends  upon 
both  the  impulse  and  the  impulse  per  unit  area,  the  later  being  the  integral 
over-time  of  the  peak  pressure.  In  the  simplest  of  cases  this  pressure  is 
determined  by  LSD  shock  wave  jump  conditions  during  the  laser  pulse  followed 
by  an  expansion  first  in  one  and  then  in  three-dimensions. 

Measurements  of  this  pressure  during  the  present  test  series  were  com¬ 
plicated  by  the  effects  of  a  finite  ignition  time  and  a  time  varying  inten¬ 
sity.  The  issue  was  further  complicated  by  the  high  intensities  used  which 
produced  shock  waves  in  which  radiation  was  important,  and  this  also  affected 
the  pressure  and  its  time  history.  The  details  of  this  early  time  history 
determine  the  impulse  per  unit  area  and  are  thus  important  for  damage. 

The  PSI  code  LSDNS  calculates  the  quasi-one-dimensional  non-steady 
flow  with  laser  heat  addition  and  real  gas  effects.  It  does  not  include  radi¬ 
ation  transport  and  therefore  cannot  make  the  transition  to  radiation-driven 
waves  expected  at  the  highest  intensities.  It  has  been  used  to  determine  the 
effects  of  the  finite  ignition  time  and  the  time  varying  intensity  on  the 
pressure  measureed  at  the  wall. 

The  time  history  of  the  laser  pulse  is  shown  in  Fig.  2.1.  The  pres¬ 
sure  behind  the  shock  wave  Ps  and  the  pressure  at  the  wall  Pw  are  shown  in 
Fig.  2.2.  The  difference  between  these  two  pressures,  a  factor  of  3,  is  in 
good  agreement  with  what  would  be  expected  from  simple  theory;  (Pw/Ps  = 
(y+1/2Y)2Y/Y-1  =  .35  where  Y  =  1.2),  and  indicates  that  the  variation  of  in¬ 
tensity  with  time  has  little  effect  on  the  peak  pressure  at  the  wall.  However 
the  measured  pressures  are  somewhat  lower  and  might  be  indicative  of  the 
importance  of  radiation  on  the  wave. 

The  effect  of  radiation  is  to  produce  higher  wave  velocities  than 
would  be  predicted  using  a  LSD  wave  criteria  and  this  is  discussed  in  more 
length  in  the  section  on  Supersonic  Radiation  Waves.  The  effect  of  a  higher 
wave  velocity  is  a  lower  pressure  behind  the  wave  and  this  is  illustrated  in 
Fig.  2.3.  It  is  important  to  note  that  the  fall  off  in  pressure  behind  the 
is  to  some  extent  compensated  for  by  the  decrease  in  the  pressure  ratio  be- 
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CHANGE  IN  LSD  WAVE  PRESSURE  WITH  VELOCITY 
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Fig.  2-3  Change  in  LSD  wave  pressure  with  velocity. 
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tween  the  wave  and  the  wall  for  small  increases  in  wave  velocity.  This  is 
because  the  faster  wave  velocity  produces  a  lower  gas  velocity  behind  the 
wave.  The  onset  of  radiation  effects  on  wave  velocity  is  best  detected  by 
measurements  of  the  wave  velocity. 

However  the  presence  of  a  lower  gas  velocity  shortens  the  communica¬ 
tion  time  between  the  shock  wave  and  the  wall.  The  pressure  increase  that 
occurs  at  the  shock  front  due  to  an  increase  in  laser  intensity,  for  example 
that  caused  by  the  secondary  pulse  occurring  at  about  0.5  ys  in  Fig.  2.1,  must 
propagate  upstream  against  the  gas  flow  to  reach  the  wall.  The  shortening  in 
time  is  evident  in  the  measurements  where  the  measured  time  of  the  pressure 
increase  at  the  wall  due  to  the  secondary  pulse  at  1  ys  (see  Fig.  2.4)  is  con¬ 
siderably  shorter  than  the  calculated  time  1.5  ys  as  shown  in  Fig.  2.2.  This 
should  be  considered  indicative  of  high  intensity  radiation  wave  phenomena, 
and  is  important  for  calculating  the  impulse  per  unit  area. 
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A  description  of  LSDNS  and  the  modifications  made  to  it  to  do  these 
calculations  is  presented  in  the  next  section  and  the  manner  in  which  we 
handled  the  real  gas  is  described  in  the  following  section. 

2 . 1  Calculation  of  Wall  pressure 

PSI  has  a  computer  code  called  LSDNS  which  calculates  quasi-one-dimen- 
sional  non-steady  flow  with  laser  heat  addition.  This  code  has  been  used  to 
calculate  the  fluid  mechanics  phenomena  occurring  in  pulsed  laser  propulsion 
thrusters,  for  example.  It  has  the  ability  to  capture  a  laser-supported  deto¬ 
nation  wave  caused  by  the  absorption  of  laser  energy.  The  thermodynamic  prop¬ 
erties  of  the  gas  which  is  absorbing  the  laser  energy  are  input  as  a  table,  so 
the  code  can  handle  gases  in  thermochemical  equilibrium.  The  absorption  coef¬ 
ficient  of  the  gas  at  laser  wavelengths  is  specified  as  a  function  of  the 
thermodynamic  state  of  the  gas. 

The  development  and  application  of  this  code  to  the  fluid  mechanics  of 
pulsed  laser  propulsion  is  described  in  Refs.  2  and  3.  In  the  present  work, 
the  code  was  used  to  calculate  the  pressure  at  x  =  0,  where  a  reflection 
boundary  condition  is  imposed.  This  boundary  condition  represents  the  condi¬ 
tion  at  a  solid  wall,  so  p  at  x  =  0  is  the  wall  pressure. 
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In  the  present  case,  the  code  was  used  in  its  purely  one-dimensional 
form.  Hie  equations  of  conservation  of  mass,  momentum  and  energy  are: 

3p  .  3pu 
3t  +  3x  * 


3pu  3 _ 

3t  +  3x 


(P 


+  pu2) 
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♦r)l 


3x 


[pu(( 


Here  p,  p,  u,  e  are  density,  pressure,  velocity  and  internal  energy  per  unit 
mass,  respectively.  I  is  the  laser  intensity,  which  obeys  the  radiation 
transport  equation  in  the  absence  of  emission. 

dl/dx  =  k^I  , 

where  kL  is  the  laser  energy  absorption  coefficient.  The  solution  of  this 
transport  equation  is 


00 

I  =  I  (t)exp (-/  k  dx)  , 

O  Li 

X 


where  the  incoming  laser  intensity  IQ  is  at  x  =  '».  Notice  that  the  sign  of 
the  transport  equation  is  chosen  so  that  the  laser  intensity  decreases  with 
decreasing  x,  since  the  laser  beam  is  propagating  in  the  negative  x  direc¬ 
tion,  toward  the  wall. 


To  complete  the  formulation  of  the  problem,  we  must  specify  the  time 
dependence  of  IQ,  the  dependence  of  kL  on  the  gas  state,  and  the  thermal  and 
caloric  equations  of  state,  p(p,  e),  T(p,  e).  Then  solution  of  the  equations 
will  provide  the  desired  wall  pressure  p( x  =  0,  t) . 


2.2  Laser  Absorption  Coefficient 

The  choice  of  the  absorption  coefficient  was  made  with  a  view  toward 


simplicity.  The  only  real  purpose  of  kL  is  to  get  the  energy  into  the  gas,  so 
an  elaborate  formulation  was  neither  necessary  or  possible  within  the  con¬ 
straints  of  the  present  effort.  A  further  reason  for  simplicity  was  suggested 


by  the  size  of  the  equation  of  state  table  that  would  be  needed.  For  the  sol¬ 
ution  of  the  flow  equations  only  p  and  T  as  functions  of  p  and  e  are  needed. 

If  kL  requires  species  concentrations ,  then  they  would  also  have  to  be  in  the 
table . 

Based  on  these  factors,  k^  was  chosen  as 
kL  =  10“3  for  T  <  5000  K 

3  pTo  2 

=  1  o  f — — )  for  T  >  5000  K  , 
p  T 
o 

where  pQ  and  T0  are  the  conditions  before  the  laser  pulse.  The  quadratic  de¬ 
pendence  on  p/T  is  suggested  by  the  dependence  of  inverse  bremsstrahlung  ab¬ 
sorption  on  the  product  of  the  electron  and  ion  number  densities. 

2.3  Equation  of  State 

Air  is  the  gas  of  interest  in  the  present  work.  The  natural  variables 
for  the  flow  equations  are  seen  to  be  p  and  e.  From  these  two  thermodynamic 
variables,  all  the  others  can  be  obtained  under  conditions  of  thermochemical 
equilibrium.  Of  course,  e  cannot  be  used  as  an  independent  variable  directly. 
The  equations  for  equilibrium  air  can  be  formulated  using  p  and  T  as  inde¬ 
pendent  variables.  Then  an  iteration  scheme  can  be  used  to  find  T  for  a  given 
p  and  e.  With  this  T,  p  and  any  other  desired  variables  can  be  found.  The 
method  of  solution  of  the  flow  equations  requires  use  of  the  local  speed  of 
sound  to  determine  the  time  step,  so  the  speed  of  sound  is  included  in  the 
equilibrium  calculation. 

During  the  solution  of  the  flow  equations,  the  expressions  T(p,  e)  and 
p(p,  e)  must  be  evaluated  many  times.  To  avoid  having  to  iterate  on  e  each 
time,  a  table  of  T,  p,  and  a  (speed  of  sound)  is  constructed,  using  p  and  e  as 
independent  variables.  This  table  is  constructed  by  the  iteration  process, 
but  the  number  of  entries  is  several  orders  of  magnitude  less  than  the  number 
of  evaluations  of  T,  p,  and  a  needed  in  the  course  of  solution.  For  these 
evaluations,  quadratic  interpolation  in  the  table  is  used. 


The • formulation  of  a  model  of  equilibrium  air  is  described  in  the  next 
section.  The  air  is  taken  as  a  mixture  of  oxygen  and  nitrogen,  with  no  trace 


elements  included.  Up  to  20,000  K,  species  included  are  N,  N+,  NO,  N0+,  N2, 
N2+,  0,  0+,  0",  02#  02+z  C>2-'  anci  E*  An  equilibrium  program  is  used  which 
utilizes  fits  to  the  partition  functions  good  up  to  20,000  K,  and  takes  p  and 
T  as  independent  variables.  Above  20,000  K,  a  different  equilibrium  program 
is  used,  which  deals  only  with  atoms  and  ions,  but  no  molecules.  It  handles  a 
mixture  of  two  ionizing  monatomic  gases,  and  in  the  present  calculations  used 
6  ions  of  each  atom,  oxygen,  and  nitrogen.  It  includes  lowering  of  the  ion¬ 
ization  potential.  It  uses  as  input  the  energy  levels  and  weights  of  the  in¬ 
ternal  state  of  each  atom  and  ion,  suitably  grouped. 

Using  this  model  of  air,  the  table  of  T,  p,  and  a  as  a  function  of  p 
and  e  was  constructed  for  use  in  the  solution  of  the  flow  equations. 

2.4  Laser  Intensity 

Laser  intensity  was  modeled  by  fitting  the  experimental  intensity  pro¬ 
vided.  A  table  was  available  providing  y/yp  vs.  time,  from  21  to  729  ns.  A 

log-log  plot  showed  that  the  curve  could  be  represented  by  five  straight  lane 
segments.  Least  square  fits  to  these  five  segments  provided  the  following 
equations: 

64  <  t  <  94: 

94  <  t  <  125: 

125  <  t  <  409 

409  <  t  <  521 

521  <  t  <  625 

Here  t  is  in  nanoseconds. 

The  fit  was  made  between  64  and  625  ns  because  the  values  after  625  ns 
appeared  to  be  zero,  and  we  allowed  the  first  70  ns  for  the  ignition  process. 
The  flow  field  code  LSDNS  was  started  at  the  70  ns  point,  and  the  laser  was 
turned  off  at  625  ns. 


y/yP  =  2.214  x  10-4  t1-8536  (1) 
y/yp  =  0.949  (2) 

y/yp  =  3.847  x  1 o5  t~2>67  (3) 
y/yp  =  4.6125  X  10“28  t9*889  (4) 
y/yp  =  1.8075  x  1020  t-7*65  .  (5) 


Equations  (1-5)  were  used  over  intervals  determined  by  the  fits,  with 
the  joining  points  chosen  where  the  fitting  equations  were  equal.  The  range 
of  t  for  each  of  the  equations  is  the  following: 


70.  < 

t  <  91 .09. 

(1  ) 

91  .09 

<  t  <  126. 

(2) 

126.  < 

t  <  422. 

(3) 

422.  < 

t  <  521  . 

(4) 

521.  < 

t  <  625. 

(5) 

In  addition  to  the  mean  values  just  modeled,  an  oscillation  was  impos¬ 
ed,  of  the  form 

F  =  1  +  Ae_Bt  cos  D  t  . 


Where  t  is  specified  in  ns.  Ifie  constant  D  was  chosen  so  as  to  produce  a  com¬ 
plete  oscillation  in  17  ns,  so  that  D  =  2tt/17  =  0.3696.  The  other  two  con¬ 
stants  were  chosen  by  specifying  that  the  amplitude  of  the  oscillation  at  the 
end  of  4  cycles  (68  ns)  was  1.79,  and  at  the  end  of  8  cycles  (136  ns),  it  was 
1.32.  These  result  in  the  values  B  =  0.01329  and  A  =  1.950.  Therefore,  F  was 
taken  as 


F  =  1  +  1.950  e“°*07329tcos  0.3696t  . 
The  final  intensity  was  then 
IQ(t)  =  KF(y/yp) 

starting  at  70  ns,  where  K  is  a  scaling  constant. 


A  plot  of  y/yp  and  IQ(t)  is  shown  in  Fig.  2.1,  with  K  =  1. 


The  total  fluence  incident  up  to  any  time  is  just  the  time  integral  of 
[Q.  Although  the  expression  can  be  integrated  analytically,  it  is  simpler  to 


just  evaluate  IQ  at  small  time  intervals  and  perform  the  summation.  The  scal¬ 
ing  constant  K  is  used  to  specify  the  fluence  desired  by  using  the  integral  of 


IQ  between  the  desired  times. 
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This  appendix  describes  the  model  used  for  high  temperature  air.  The 


model  involves  two  different  temperature  regimes,  above  and  below  20,000  K. 

Below  20,000  K  PSI  has  a  standard  equilibrium  program  which  uses  curve 
fits  to  the  partition  functions  to  determine  the  equilibrium  constants  for  a 
gas  mixture.  For  use  below  6000  K,  the  data  in  the  JANAF  tables  are  used. 

For  certain  species,  curve  fits  applicable  up  to  20,000  K  have  been  construct¬ 
ed.  These  fits  were  made  in  a  form  compatible  with  the  analytical  expressions 
for  the  partition  functions.  The  natural  logarithm  of  the  partition  function 
Q  was  expressed  as 


»  »  i  i  2  » 

£nQ  =  L1  +  L^lnT  -  L3/T  +  I^/T  +  LgT 
'2  '3 

+  L  T  +  L_T  +  v(T) 
b  / 


The  constants  were  determined  by  a  combination  of  analytical  expres¬ 
sions  and  curve  fitting.  At  moderate  temperatures,  Q  is  a  combination  of  ro¬ 
tational,  vibrational  and  electronic  partition  functions: 


Q  -  QtQrQvQe • 


The  expressions  for  the  component  partition  functions  are: 


Qt  =  (2rmkT/hp) 3/2 
Qr  =  T/o9r 


Qv  -  n 


e-0v/T)-1 


Qe  =  ge-°/T 


Here  m  is  the  species  mass,  k  is  Boltzmann's  constant,  hp  is  Planck’s  con¬ 
stant,  9r  is  the  rotational  temperature  for  diatomics,  o  is  unity  except  for 
homonuclear  molecules  for  which  a  =  2,  9V  is  the  vibrational  temperature  for 
diatomics,  g  is  the  degeneracy  of  the  ground  electronic  state,  and  D  is  the 
energy  of  the  species  at  zero  temperature,  expressed  in  K. 


The  constants  Ig  -  L3  and  v(T)  were  used  to  agree  with  these  simple 


expressions: 


Ig  “  1.5  An(2nmk/hp)  +  £n(g/a0r) 


=  2.5  for  diatomic  species 


1.5  for  monatomic  species 


L3  =  D 

v(T)  =  -£n(1  -  e“®v/T)  for  diatomic  species 
=0  for  monatomic  species 

I  I 

The  remaining  constants  L4  -  L7  were  used  to  fit  the  expression  for  £nQ  to  a 
calculation  of  Q  which  included  many  terms  in  the  electronic  partition  func¬ 
tion,  anharmonicity  in  vibration,  and  vibration-rotation  coupling,  but  not 
lowering  of  the  ionization  potential.  This  method  of  fitting  provided  good 

I  ( 

fits,  and  also  small  values  of  the  constants  L4  -  L7. 

Air  is  taken  to  be  a  mixture  of  oxygen  and  nitrogen,  in  the  molar 
composition  of  78.849%  N2,  21.151%  02,  a  ratio  of  0.26825  02  to  each  N2.  The 
species  which  have  been  fitted  include  N,  N+,  NO,  N0+,  N2,  N2+,  0,  0+,  0",  02, 
02+ ,  02",  E. 


Above  20,000  K,  there  are  no  molecules  present.  We  assume  that  only 
atomic  species  and  their  positive  ions  are  present.  Therefore,  air  is  repre¬ 
sented  as  a  mixture  of  2  monatomic  ionizing  gases.  A  computer  program  for  the 
equilibrium  state  of  such  a  mixture  was  developed,  based  on  a  previous  program 
for  a  single  monatomic  ionizing  gas. 


'•'iTiliV*? 


We  consider  a  mixture  of  two  gases,  denoted  by  superscripts  1  and  2. 

For  each  gas,  there  are  atoms  and  N-1  ions,  plus  electrons,  a  total  of  2N  +  1 

1  12  2  12 

species.  The  number  densities  are  nA  =  n^ ,  nA  =  n^ ,  n^,  n^(i  =  2,3,..N)  and 

1  12  2 

nE.  The  corresponding  masses  are  mA  =  m^  ,  mA  =  m-|  ,  =  mA  -(i  -  1)mE  for 

both  species  1  and  2,  and  mE.  The  charges  are  0,  (i  -  1)  and  -1  electronic 
charges . 


The  total  number  density  is 


12  N  i  2 

n  =  nA  +  nA  =  nE  +  i|2  +  nj 


charge  conservation  requires 


N 


-n  +  £  (i  -  1  )  (n^  +  n2) 

E  i=1  11 


The  density  of  the  mixture  is 

11  2  2 


11  i  £  „  r  1  .  1  .  . 

p  =  ran.  +  inn  +  m  n„  +  £  n.(m  -  (i 

A  A  AA  EE  „  L  l  A 


i=2 


1  )m  ) 
E 


+  ni(mA  "  (i  "  1)mE)J 


The  use  of  charge  conservation  simplifies  this  to 

N  „  _  N 


1,1  r  1,  2,2  _  2, 

p=m^(n.  +  £  n,l  +  m, (n  +  £  n.) 

AA  .  „  l  AA  .  „  l 

1=2  1=2 


The  mole  ratio  of  the  mixture  is  defined  to  be 


2  N  2  1  N  1 

f  =  (n,  +  £  n; )/(n'  +  £  n!) 
n  A  .  _  l  A  .  _  l 

1=2  1=2 


(2.1  ) 


(2.2) 


(2.3) 


(2.4) 


which  allows  the  total  number  density  to  be  written  as 


n=(nl+£  n] ) ( 1  +  f  )  +  n 


(2.5) 


If  we  also  define  the  mass  ratio  as 


f  =  nu/m][  (2.6) 

m  A  A 

the  density  Eq.  (2.3)  becomes 

p  =  mhn  -  n  )  <1  +  f  f  )/( 1  +  f  )  (2.7) 

A  E  n  m  n 


The  gas  law  is 


P 


nkT  = 


pT 


m ' 


-  Vn 


1  ~  fn 

+  f  f 
n  m 


(2.8) 


The  equilibrium  state  is  specified  by  the  Law  of  Mass  Action  relating 
the  partial  pressures  of  the  i-th  and  (i+1)-th  ions  to  the  electron  partial 
pressure: 


p  p.  + 


VP, 


=  f.  (T) 
IP 


Here  f^p  is  the  equilibrium  constant  for  the  ionization  of  the  i-th  ion,  which 
will  be  defined  later. 


If  we  introduce  the  normalized  reaction  variables 

a.  =  n./n  =  P./P 
x  1  1 

then  Mass  Action  becomes 


a  a.  ./a- 
E  x+r  l 


f . 
ip 


/P 


(2.9) 


(2.10) 


LUC  UOC 


which  is  the  form  appropriate 
For  the  present  application,  we  need  p,  T  as  independent  variables,  so  we  re¬ 
write  Eq.  (2.10)  using  the  gas  law  Eq.  (2.8): 

1 


1  -  a 


* .  ,  f .  m '  1  +  f  f 
1+1 _ i£  _A  _ n  m 

a.  PT  k  1  +  f 


(2.11) 


n 


There  are  N-1  ions  of  each  species,  and  an  equation  like  this  for  the 
formation  of  each  ion.  Thus  there  are  2 ( N— 1 )  such  equations. 
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..... 


In  terms  of  the  a  variables,  the  number  density  relation,  Eq.  (2.1), 
the  change  conservation  relation,  Eq.  (2.2),  and  the  mole  ratio  relation,  Eq. 
(2.4),  may  be  written 


.  1  2  1.12, 
1  =  a  +  a  +  a  +  I  (a.  +  a.) 
A  A  E  ,  _  l  l 
1=2 


-a  +  E  (i  -  1 ) (a!  +  a?)  =  0 
E  i=2  11 


2  *2,1  ? 

a  +  E  a.=(a  +  E  a.)f 

A  .  „  i  A  .  _  l  n 

1=2  i=2 


(2.12) 


(2.13) 


(2.14) 


The  lowering  of  the  ionization  potential  is  calculated  using  the 
Debye-Huckel  theory.  The  decrement  in  the  potential  is  given  by 


_  1 /2  N  -  ~  _  1/9 

A  I.  =  2ie  (— )  (  E  (n.  +  n , )  ( j  —  1 )  +  n  )  ' 

1  kT  j=2  3  3  E 


(2.18) 


where  e  is  the  electron  charge  in  esu.  If  Al^  is  expressed  in  ergs,  2eVu/k  = 
3.343  x  10~20,  while  if  it  is  expressed  in  cm'1,  one  divides  by  hpc  so  the 
constant  is  2e3Ar/k/hpc  =  1  .683  x  10~4. 

1  2 

Since  Al^  involves  n^,  n^  and  ng,  while  Ig^  is  needed  on  the  right 
side  of  Eq.  (2.18)  to  calculate  these  quantities,  the  solution  must  clearly  be 
iterative.  However,  solving  for  the  a  is  already  iterative.  The  procedure 
followed  was  simply  to  use  the  current  values  of  the  number  densities  to  cal¬ 
culate  the  Al^  used  to  find  the  next  iteration  on  a.  In  view  of  the  approxi¬ 
mate  nature  of  the  Al^  correction,  this  procedure  is  sufficiently  accurate. 

The  internal  energy  per  unit  mass  of  each  species  is  expressed  in 


terms  of  the  internal  partition  function  by 


kT  r  3 

en  =  —  —  + 


•  i  3  kT 

‘J '  eE  =  o  =  — 


(2.19) 


where  the  first  term  is  the  translational  part.  The  mixture  internal  energy 
is  then 


e  =  E  e^/p  +  egPg/P/ 


(2.20) 


Using  Bq.  (2.16)  in  Eq.  (2.19)  and  performing  the  sum  in  Eq.  (2.20),  we  find 


The  sum  over  i  is  carried  out  over  all  the  atoms  and  ions,  2N  terras,  but  not 
over  the  electrons  which  have  no  internal  energy. 


The  enthalpy  per  unit  mass  is  found  from  e  as 
h  =  e  +  p/p 


(2.22) 


The  formulation  so  far  has  treated  multiple-ionizing  atoms.  If  we  are 
interested  in  using  this  model  to  represent  the  high  temperature  state  of  a 
mixture  of  diatomic  molecules,  where  there  are  no  molecules  present,  the  ques¬ 
tion  of  the  energy  reference  level  arises.  So  far,  the  energy  reference  level 
has  been  that  of  the  atoms  at  zero  degrees  K.  However,  in  most  application 
to  diatomic  molecules  the  reference  level  used  is  the  molecules  at  zero  de¬ 
grees  K.  This  requires  the  addition  of  energy  to  the  expressions  for  e  and  h, 
to  refer  them  to  the  correct  level.  This  is  simply  done  by  adding  to  each 
atom  and  ion  one  half  the  dissociation  energy  per  molecule.  If  these  disso¬ 
ciation  energies  are  denoted  by  D1 ,  D2,  then  the  internal  energy  and  enthalpy 
referred  to  the  molecule  is  related  to  the  values  in  Eqs.  (2.21)  and  (2.22) 


a  1  A  A 

e*  =  e  +  E  (n.D  +  n^D  )/2p 
,  l  l 

i  =  1 


(2.23) 


h*  =  h  +  E  (nV  +  n2D2)/2p 
i=1  1  1 


(2.24) 


E  g. .e . .e"£ij/kT 

3  P  +  e  iH  [£*  +  0— - ] 

2  p  .  P  L  i  z  e-eijAT  J 

j  ij 


(2.21 ) 


The  speed  of  sound  can  be  found  from  the  specific  heats  and  deriva¬ 
tives  of  p  with  respect  to  the  independent  variables  p  and  T.  The  specific 
heat  at  constant  volume  is 


cv  =  (3e/3T)p 

while  the  specific  heat  at  constant  pressure  is  related  to  cv  by 


2.25) 


'P  -  p 
The  speed  of  sound  squared  is 


CP  -  cv  +  77  (3T}p/l3p  t 


(2.26) 


|2  .  !e  ,i£ 


( 3  )_ 

V  dp  T 


(2.27) 


These  relations  express  cp,  cv  and  a2  entirely  in  terms  of  derivatives  which 
use  p  and  T  as  independent  variables.  These  derivatives  can  be  found  by 
finite  difference  expressions  once  e  (or  e* )  and  p  are  known  from  Eqs .  (2.21) 
(or  2.23)  and  (2.3).  The  same  finite  difference  representations  of  the  deriv¬ 
atives  can  be  used  with  the  low  temperature  model  to  find  the  specific  heats 
and  sound  speeds. 


With  these  models,  the  equilibrium  state  of  air  (the  mixture  of  oxygen 
and  nitrogen)  can  be  obtained  as  a  function  of  p,T.  In  particular,  p,  e  and  a 
can  be  found,  as  well  as  the  number  densities  of  the  atoms  and  all  ions  one 
chooses  to  include. 


To  obtain  the  state  as  a  function  of  p  and  e  (or  e*),  one  can  perform 
an  iteration  on  T  until  the  desired  value  of  e  or  e*  is  found.  This  was  the 
procedure  used  to  construct  a  table  of  p,  T  and  a  using  p  and  e*  as  indepen¬ 
dent  variables.  This  table  was  used  to  supply  the  equation  of  state  for  the 
LSDNS  Code. 
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3.  IMPULSE  SCALING 


Although  the  pressure  time  history  on  the  target  during  and 
immediately  after  the  laser  pulse  is  very  much  dependent  on  the  details  of  the 
laser  interaction,  LSD  or  radiation  wave  for  example,  the  pressure  at  late 
times  depends  mainly  on  the  total  energy  absorbed.  Thus  the  impulse,  the 
integral  over  time  of  the  force,  becomes  independent  of  the  laser  interaction 
details  in  many  cases,  and  valuable  insight  into  the  scaling  with  density  and 
energy  of  the  impulse  can  be  obtained  from  quite  simple  ideas. 

It  will  be  found  that  these  simple  ideas  lead  to  the  scaling  laws  for 
the  impulse  as  measured  by  the  pendulum.  It  is  important  to  note  here  that 
Lindbergs1  work  on  the  damage  thresholds  for  re-entry  vehicles  found  that  the 
impulse  was  not  the  only  determinant  of  failure  but  that  the  impulse  per  unit 
area  was  also  important.  He  found  that  if  the  level  of  impulse  required  for  a 
given  spot  area  was  determined  and  then  the  spot  area  was  reduced  by  say  a 
factor  of  ten  the  impulse  required  did  not  stay  the  same.  It  did  not  go  down 
by  ten  as  might  be  expected  if  it  was  mainly  the  impulse  per  unit  area,  the 
integral  of  pressure  over  time,  that  was  important;  rather,  it  went  down  a 
factor  of  three.  They  are  both  equally  important.  The  insight  these  damage 
experiments  give  us  about  what  to  expect  for  impulse  per  unit  area  measure¬ 
ments  will  be  discussed  later  when  we  discuss  scaling. 

Consider  the  geometry  shown  in  Pig.  3.1  consisting  of  a  target  disk  on 
which  the  impulse  is  to  be  determined,  imbedded  in  an  infinite  plane.  Assume 
that  the  laser  spot  size  is  small  compared  to  the  disk  radius,  and  also  that 
the  distance  the  LSD  or  radiation  wave  travels  in  the  laser  pulse  time,  ctp, 
is  also  small  compared  to  the  disk  radius.  We  approximate  the  laser  interac¬ 
tion  by  a  spherical  point  detonation  in  the  center  of  the  disk  in  the  half 
plane . 


If  the  energy  ahead  of  the  shock  wave  can  be  neglected,  Sedov2  has 
shown  that  the  solutions  for  the  flow  are  self  similar,  meaning  that  the  pro¬ 
files  of  p,  p  and  T  behind  the  wave  change  in  scale  as  the  wave  moves  but  not 
in  shape;  that  is  the  pressure  p  is  such  that  P/ps  =  f(r/rs)  where  ps  the 
pressure  at  the  shock  point  and  the  rg  position  of  the  shock  point  vary  with 


5773-6 


1  ATM  AIR 


SHOCK  WAVE 


TARGET 


Contribution  to  impulse  comes  from  time  when  shock  wave  is 
on  target  plus  time  when  it  is  off  target 


Fig.  3.1  Laser  intensity  used  in  wall  pressure  calculation. 


time  but  the  function  f  does  not  change.  As  the  shock  wave  continues  to  ex¬ 
pand  the  energy  per  unit  mass  in  the  gas  behind  the  wave  decreases  and  a  point 
is  reached  where  the  energy  per  unit  mass  of  the  gas  ahead  of  the  shock  is  no 
longer  negligible  and  the  spatial  profiles  will  depend  upon  the  energy  of  the 
gas  ahead  of  the  shock  wave.  The  distance  at  which  this  occurs  can  be  esti¬ 
mated  by  setting  the  energy  in  the  original  point  explosion  equal  to  the 
energy  in  the  background  gas  the  shockwave  has  engulfed. 


4irr' 


P1  =  Eo 


and  thus  we  define  a  distance  rQ  by 


E  1/3  r 

o  .  s 

r  =  —  and  Z  =  — 

o  P,  r 

1  o 


such  that  the  effects  of  the  energy  ahead  are  important  when  Z  is  of  order  1 . 
The  change  in  the  self-similar  pressure  profile  for  Z  near  one  has  been  given 
by  Sedov2  and  is  reproduced  in  Fig.  3.2.  An  actual  pressure  curve  using  the 
curves  and  table  of  Fig.  3.2  is  shown  in  Fig.  3.3  for  the  value  of  Z  =  .7566. 
Notice  that  the  pressure  at  the  origin  is  negative,  that  is  below  atmospheric. 
At  this  time  the  impulse  on  the  target  is  decreasing  with  time.  This  will  be 
shown  to  be  an  important  effect  for  large  targets,  that  is,  targets  whose  size 
is  comparable  to  rQ  the  energy  radius.  This  effect  is  probably  not  important 
for  cases  of  interest  for  damage  because  the  pressures  are  quite  low,  approx¬ 
imately  the  ambient  pressure.  They  are  important  for  interpreting  experi¬ 
mental  data;  for  example  all  of  the  10.6  y  single  pulse  data  was  in  this 
regime . 


For  large  energies  simple  blast  wave  scaling  is  appropriate.  This 
scaling  conserves  total  energy.  Combining  this  with  conservation  of  momentum 
across  the  shock  wave  the  impulse  scaling  laws  can  be  obtained. 


Conservation  of  energy  pr3  «,  E 
Conservation  of  momentum  pc2  a  p 


where  p  is  the  pressure  behind  the  wave 
r  is  the  position  of  the  wave 
E  is  the  energy  put  in  by  the  laser 
p  is  the  gas  density  ahead  of  the  wave 
c  is  the  velocity  of  the  wave 

To  obtain  the  impulse  on  the  target  we  wish  to  evaluate 


I  =  f  pAdt 


In  this  integral,  A  is  the  target  area,  and  the  integral  is  best  evaluated  by 
considering  two  times,  early  times  when  the  shock  wave  is  still  on  the  target 
and  late  times  when  it  is  off  the  target.  To  do  this  it  is  easiest  if  the 
subs  ti tution 


dr 

dt  =  — 
c 


is  made, 


I  =  /  PA 


r  C 

J  PA- 


where  rT  is  the  target  radius 


T 
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2  dr  r  , E  , 
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2  dr 
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pr 


Notice  both  integrals  are  such  that  their  major  contribution  derives  from  the 
time  the  blast  wave  radius  is  equal  to  the  target  radius.  One  can  then 
write, 

3/2  1/2  1/2 
I  «  ,93  r  E  p 

T 

where  the  constant  has  been  derived  from  more  detailed  calculations 2 > 4  .  Thus 
we  have,  that  as  long  as  we  stay  with  large  laser  energies,  the  impulse  cou¬ 
pling  coefficient  I/E  is  expected  to  decrease  as  the  inverse  square  root  of 
the  energy  and  as  the  square  root  of  the  density. 


If  an  attempt  is  made  to  evaluate  the  impulse  per  unit  area  in  a  simi¬ 
lar  manner  the  integral  diverges,  the  high  blast  wave  pressure  at  early  time 
dominates.  In  contrast  to  the  case  for  the  impulse  where  the  early  time  de¬ 
tails  were  unimportant,  the  early  time  details  determine  the  impulse  per  unit 
area.  Remember  damage  depends  on  both  of  these  quantities.  PSI  has  examined 
the  blast  wave  scaling  in  more  detail  so  as  to  avoid  this  difficulty  and  at 
the  same  time  remove  the  earlier  detailed  restrictions,  rs  <  rT  and  CTp  <  rT. 
This  work  is  described  in  the  next  section.  This  enables  us  to  determine  the 
impulse  under  more  general  conditions  and  also  to  determine  the  impulse  per 
unit  area  and  how  it  scales.  However  the  PSI  theory  restricts  itself  to  the 
strong  shock  limit;  negative  pressures  are  not  considered  as  they  are  not 
interesting  for  damage  as  mentioned  earlier.  Koslova,  et  al.3  have  carried 
the  simple  spherical  blast  wave  analysis  into  this  region  and  these  results 
can  be  used  to  interpret  the  10.6  pm  data  which  was  taken  in  this  regime. 

They  found  that  using  an  ideal  gas  and  a  point  explosion  recreates  the  general 
effects  but  were  higher  than  the  data  by  a  factor  of  two.  They  were  required 
to  use  a  distributed  explosion  and  real  gas  effects  to  reproduce  the  data. 

An  interesting  question  is  how  much  less  than  rT,  rs  has  to  be  in 
order  to  ignore  it  and  approximate  the  solution  by  a  spherical  blast  wave 
solution.  For  the  impulse  per  unit  area  the  answer  is  never,  as  was  stated 
before.  For  the  impulse  the  data  and  PSI  theory  indicate  that  even  for  laser 
spot  sizes  equal  to  the  target  radius  the  simple  spherical  blast  wave  solution 
is  good  enough.  For  ct  we  have  no  data  as  cr  was  always  much  less  than  the 
target  radius. 

Figure  3.4  shows  the  1 .06  p  data  taken  at  Battelle  on  aluminum  and 
compares  it  to  the  PSI  theory  appropriate  for  high  laser  energies,  to  the  the¬ 
ory  of  Kozlova  and  also  to  the  Russian  data  which  has  been  summarized  by  the 
dashed  curve.  Note  that  the  effect  of  the  late  time  negative  processes  is 
noticeable  at  values  of  rT/rQ  as  low  as  .2.  Figure  3.5  shows  the  data  for 
Grafoil.  Figure  3.6  compares  the  10.6  p  impulse  data  with  the  Russian  1.06  p 
data.  The  laser  energy  has  been  multiplied  by  0.56  for  the  10.6  p  data  to 
account  for  our  best  estimate  as  to  the  absorbed  portion.  Notice  that  this 
data  is  predominately  in  the  large  target  regime. 
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4.  SUPERSONIC  RADIATION  WAVES 


When  a  piston  travels  down  a  cylinder  filled  with  a  gas,  it  pushes  the 
gas  ahead  of  it  compressing  and  heating  it.  The  compression  and  heating  is 
accomplished  by  gas  collisions  in  a  narrow  region  known  as  a  shock  wave.  If 
the  piston  travels  at  a  constant  velocity,  the  shock  wave  moves  ahead  of  it  at 
a  constant,  higher  velocity  creating  a  region  of  compressed  heated  flow.  The 
properties  of  the  gas  in  this  region  can  be  determined  by  conserving  mass,  mo¬ 
mentum  and  energy  across  the  wave.  For  strong  shocks  where  the  gas  energy 
ahead  of  the  wave  can  be  neglected  the  conservation  equations  are: 

u  =  (l  -  )  Mass  (4.1) 

p  =  p0Du  Momentum  (4.2) 

u2 

P0D  (e  -  —  )  =  JQ  Energy  (4.3) 


e 


-  1  P 

e(p,p)  =  y-1  p 


Gas  Law 


(4.4) 


Where 


p0  =  Gas  density  ahead  of  the  shock  wave. 

D  =  Shock  wave  velocity. 

u  =  Gas  velocity  behind  the  wave  in  laboratory  coordinates, 
p,  p,  e  =  Pressure,  density  and  internal  energy  behind  the  wave. 

Y  =  The  ratio  of  specific  heats  and  is  defined  by  the  second 
equality. 

JQ  =  Represents  an  energy  input  per  unit  time  and  area  from  an 
external  source,  i.e.,  a  laser. 

For  the  case  of  the  piston,  the  gas  velocity  u  is  equal  to  the  piston 
velocity,  JQ  =  0  and  the  source  of  energy  compressing  and  heating  the  gas  is 
the  work  being  done  by  the  piston  W  =  pu.  if  the  piston's  velocity  changes  to 
a  new  constant  velocity,  this  fact  is  communicated  to  the  shock  front  at  the 
speed  of  sound  and  eventually  the  shock  strength  adjusts  to  the  new  piston 
velocity.  There  is  no  difficulty  doing  this  as  u  +  a  >  D,  that  is  sound  waves 
leaving  the  piston  at  velocity  a  can  overtake  the  shock  wave. 
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If  another  source  of  energy  (chemical  or  from  a  laser)  exists,  (JQ  *  0) 
this  simple  connection  between  the  piston  velocity  and  the  shock  velocity  is 
lost.  There  is  an  unsteady  flow  field  between  the  shock  and  the  piston  and  the 
velocity  at  the  piston  is  not  usually  the  same  as  the  velocity  u,  just  behind 
the  shock  wave.  With  u  no  longer  known  the  four  equations  now  have  five  un¬ 
knowns  and  another  relation  is  required.  The  possible  solutions  to  the  conser¬ 
vation  equations  are  shown  with  p  and  PQ/p  as  coordinates  in  Fig.  4.1  for  air 
with  an  absorbed  laser  intensity  J  of  2  x  109  w/cm2.  For  comparison  the  same 
curve  for  J  =  0  is  shown  in  the  lower  left  hand  corner.  (Notice  that  the  ver¬ 
tical  scale  has  been  divided  by  100) . 


Combining  Eq.  4.1  and  4.2  the  relation  p  =  p0D2  (l  -  is  obtained. 

Thus  the  slope  of  a  line  starting  at  (0,  1)  in  the  p,  p0/p  plane  measures  the 
quantity  p0D2.  The  minimum  value  of  D  is  given  by  the  line  shown  tangent  to 
the  curve  at  A.  Jouquet  (1905)  pointed  out  this  minimum  velocity  satisfies  the 
condition  that 


u  +  a  =  D.  (4.5) 

Strong  detonations,  corresponding  to  intersection  points  to  the  left  of 
A,  have  u  +  a  >  D,  thus  information  generated  at  the  piston  can  reach  the  shock 
front;  these  provide  solutions  to  the  flow  for  large  piston  velocities  up  >  uA 
where  u^  is  the  gas  velocity  at  the  point  A.  Weak  detonations,  corresponding 
to  intersection  points  to  the  right  of  A,  have  u  +  a  <  D  and  imply  an  indeter¬ 
minacy  in  specifying  the  flow  for  piston  velocities  up  <  u^.  (Usually  the  gas 
piston  velocity  is  zero  for  the  cases  of  interest  for  laser  target  interac¬ 
tions.)  Chapman  and  Jouquet  hypothesized  that  the  wave  would  not  run  in  this 
indeterminate  region  but  at  the  point  A  and  these  waves  satisfying  Eq.  4.5  are 
called  Chapman-Jouquet  detonation  waves.  Using  arguments  based  upon  a  combus¬ 
tion  process  propagation  mechanism,  von  Neuman  (1942)  demonstrated  that  the 
weak  detonation  waves  were  in  fact  impossible.  If  the  waves  are  to  run  in  this 
region  a  propagation  mechanism  other  than  collisions  (the  combustion  process) 
is  required. 


The  addition  of  Eq.  (4.5)  to  the  original  set  of  four  now  completes 
the  set  and  the  equations  can  be  solved.  The  gas  velocity  just  behind  the  wave 
can  be  determined  and  the  unsteady  flow  field  between  the  wave  and  the  piston 
can  be  calculated. 
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When  collisions  are  the  mechanism  of  propagation  for  the  wave,  use  of 


Eq.  4.5  closes  the  set  of  equations.  However,  if  another  propagation  mech¬ 
anism  exists  the  Hugoniot  curve  shown  in  Fig.  4.1  admits  higher  velocity 
solutions  such  as  the  one  shown  for  point  B.  Notice  that  the  faster  solutions 

have  1  and  from  Eq.  (4.1)  this  implies  that  the  gas  velocity  is  close 

P 

to  zero.  This  is  a  wave  that  goes  very  quickly  through  the  gas,  raising  its 
pressure,  density  and  velocity  only  slightly.  Radiation  can  provide  a  mech¬ 
anism  for  propagation  of  such  a  wave. 


The  wave  absorbs  the  laser  radiation,  which  is  coming  from  in  front  of 
the  wave,  producing  a  hot  gas  of  temperature  Tf  with  internal  energy  Ef.  This 
hot  gas  radiates  energy  forward  (and  backward),  heating  the  gas  in  front  of  the 
wave.  The  heated  gas  becomes  ionized  and  can  in  turn  absorb  the  laser  radi¬ 
ation  through  inverse  bremsstrahlung,  thus  providing  a  mechanism  for  propagat¬ 
ing  the  wave  forward.  Although  in  front  of  the  wave,  the  unheated  gas  cannot 
absorb  the  laser  radiation,  but  does  absorb  the  radiation  from  the  hot  gas  be¬ 
hind  the  wave  which  radiates  predominantly  in  the  deep  ultraviolet.  The  phy¬ 
sics  just  outlined  can  be  expressed  by 

p0De*  =  eaTf^  (4.5A) 

where 

e*  is  the  energy  required  by  the  air  in  front  of  the  wave  before  it 
can  absorb  the  laser  radiation 


e  is  the  emissivity  of  the  laser  heater  gas 

o  is  the  Stef an-Boltzmann  constant 

Bergel’son  et.  al.1  points  out  that  as  all  the  laser  energy  is  assumed 
to  be  absorbed  at  the  wave,  there  will  be  an  overshoot  in  energy  at  the  wave 
front  equal  to  e*,  the  energy  radiated  ahead,  so  that  the  energy  right  at  the 
wave  front  ef  is  given  by 

* 

ef  =  e  +  e 

where  as  before  E  is  the  internal  energy  behind  the  wave.  Such  a  situation  is 
shown  in  Fig.  4.2.  Notice  that  we  have  used  the  fact  that  the  gas  velocity 
goes  to  zero  to  neglect  its  square  in  the  energy  equations  below.  The  energy 
Eq. (4.3)  can  now  be  rewritten 


RADIATION  FORWARD  DRIVES  WAVE 


P0D  (ef  -  e  )  =  JQ 


(4.3A) 


The  set  of  Eqs.  (4.1,  4.2,  4.3a  &  4.4),  plus  Eq.  (4.5A)  makes  a  complete  set 
and  can  again  be  solved  for  the  five  unknowns  (e^,  p,  p,  u  and  D) .  The  new 
propagation  mechanism  equation  5A  has  replaced  equation  (5).  Notice  that  for 
these  waves  u  +  a  <  D,  collisions  no  longer  can  play  a  role  in  the  propagation 
mechanism. 

The  solution  of  this  last  set  of  equations  has  been  considerably  sim¬ 
plified  by  neglecting  the  square  of  the  gas  velocity  in  the  energy  equation  in 
as  much  as  only  the  two  equations,  3A  &  5A  remain  to  be  solved  simultaneously. 
The  determination  of  the  energy  e*  depends  upon  the  absorption  properties  of 
the  heated  air,  and  thus  depends  upon  a  length  scale  Z  and  the  gas  density. 

For  1.06  p  radiation  the  absorption  is  due  to  inverse  bremsstrahlung . 

An  absorption  length  of  1mm,  a  reasonable  scale  for  waves  traveling  a  centi¬ 
meter  or  more  requires  a  gas  temperature  of  2eV  at  one  atmosphere  density. 

Using  e*  as  2eV  and  e  =  1  the  velocity  and  temperature  as  a  function  of  laser 
intensity  has  been  determined  from  equations  3A,  5A  and  the  gas  equation  of 
state.  In  practice,  since  it  is  e(T,p)  that  is  calculated  for  real  gases  and  p 
is  fixed  in  this  approximation,  a  value  of  Tf  is  chosen.  D  is  determined  from 
equation  5A  and  the  required  laser  intensity  from  3A.  The  results  of  such  a 
calculation  for  real  air  at  one  atmosphere  for  1.06  p  radiation  are  shown  in 
Figs.  4.3  and  4.4,  where  they  are  compared  to  results  for  an  LSD  wave.  These 
results  are  compared  with  velocity  data  from  the  Batelle  test  series2  in  Fig. 


Calculations  for  a  10.6  p  wavelength  are  similar.  However,  as  the 
electron  density  required  for  absorption  is  proportional  to  1 / X 2  a  a  signifi¬ 
cantly  lower  gas  temperature  is  required  at  10.6  p.  Previous  work  on  ignition3 
indicated  that  when  the  gas  temperature  reaches  0.5  eV  the  absorption 
rises  rapidly,  and  this  temperature  has  been  chosen  for  T*  for  the  10.6  p  cal¬ 
culations  . 


There  is  one  further  complication  that  must  be  addressed  before  com¬ 
paring  with  the  10.6  p  results  at  NRL.  The  NRL  laser  pulse  was  partially  mode 


locked  giving  rise  to  a  modulated  laser  intensity.  The  radiation  level  is  very 
sensitive  to  this  variation  as  it  scales  the  gas  temperature  to  the  fourth 
power.  In  evaluating  this  effect  it  is  assumed  that  the  fluctuations  are  suf¬ 
ficiently  fast  (if  ~  10  ns)  compared  to  the  propagation  time  scale,  which  is  ap¬ 
proximately  the  absorption  length  of  the  hot  gas  radiation  in  the  air  out  front 

1 0“ ^  cm 

divided  by  the  wave  velocity  ~  - Q5'cm^s '  that  the  wave  velocity  does  not  react 

to  the  fluctuations  instantaneously  but  rather  only  through  a  perceived  higher 

emissivity.  As  the  absorption  length  in  the  air  has  only  been  estimated  this 

may  be  a  poor  approximation  and  should  be  considered  as  the  maximum  effect 

these  fluctuations  can  have.  We  define  an  effective  emissivity  multiplier  as 
<Tf4> 

£m  =  <T  >4  and  evaluate  these  time  averages  by  approximating  the  energy  equation 
as  p0D  ef  =  J0  where  ef  and  JQ  are  functions  of  time  but  not  p0D.  ef  is  related 


to  Tf  by  an  approximate  equation  of  state 
ef  =  CvTf3/2 


Thus  Tf  ~  JQ2/3  and  em  = 


<J  «/3> 


If  we  assume  a  pulse  shape  as  shown  in 


Fig.  4.2b  then  em  =  2.22.  With  this  value  of  em,  a  value  of  e  of  1 ,  and  e  = 
0.5  eV  the  velocity  and  temperature  for  a  10.6  y  radiation  driven  wave  has  been 
calculated  and  is  shown  in  Figs.  4.5  and  4.6.  This  calculation  for  velocity  is 
compared  with  data  taken  on  the  NRL  device  in  Fig.  4.7.  The  analysis  of  the 
data  from  NRL  is  described  in  Section  4.1. 


The  calculations  following  Bergel'son  et.  al.  just  described  for  1.06  y 
and  10.6  y  may  be  considered  as  placing  approximate  limits  on  the  intensity  re¬ 
gion  in  which  supersonic  radiation  waves  are  expected  to  be  found,  in  as  much  as 
the  10.6  y  predictions  were  performed  for  a  minimum  value  for  T*  and  what  is 
probably  a  maximum  value  for  em  and  the  1.06  y  calculation  prediction  was  per¬ 
formed  for  what  is  probably  a  maximum  value  for  temperature.  The  measured  velo¬ 
city  for  the  1.06  um  data  is  in  good  agreement  with  theory  indicating  that  radi¬ 
ation  is  the  impor-tant  propagation  mechanism  in  this  intensity  region.  The 
10.6  y  data  again  leads  to  this  conclusion,  hut  here  the  agreement  is  not  so 
good  with  the  measured  velocities  even  fasten  than  those  predicted  using  this 
simple  theory. 


4.1 


Analysis  of  Streak  Photograph  Data 


The  analysis  of  the  streak  photograph  data  began  with  enlarging  the 
streak  photographs  as  shown  in  Fig.  4.9.  As  seen  the  luminous  front  moves  at  a 
very  rapid  velocity  at  early  times,  while  the  laser  is  on,  followed  by  a  con¬ 
siderably  lower  velocity  at  later  times,  when  the  laser  is  off.  This  is  shown 
more  clearly  in  Fig.  4.10  which  plots  distance  versus  time  for  the  first  few 
hundred  ns  from  the  photograph  of  Fig.  4.9. 

_  w 

For  two  of  the  runs  at  the  highest  intensity,  I  =  2x  1 09  — 7  ,  the 

cm * 

distance  that  the  front  moved  during  the  laser  pulse  was  sufficient  that  the 
derivative  of  the  distance  versus  time  curve  could  be  taken  and  such  a  velocity 
plot  is  shown  in  Fig.  4.11.  For  comparison  the  laser  intensity  is  also  shown  in 
this  figure.  There  is  no  time  index  on  the  streak  pictures  so  that  the  time 
axis  for  the  velocity  has  been  arbitrarily  positioned  to  place  the  maximum  vel¬ 
ocity  at  the  point  in  time  of  maximum  intensity.  This  positioning  would  lead  to 
an  ignition  time  of  40  ns,  a  very  reasonable  value. 

For  runs  at  lower  intensity  an  approximate  maximum  slope  was  measured  as 
well  as  the  position  of  the  front  at  a  fixed  time,  taken  to  be  200  ns.  For  the 
lowest  intensities  only  the  position  of  the  front  at  the  fixed  time  could  be 
determined  with  any  degree  of  accuracy.  Therefore  the  position  of  the  front  at 
fixed  time  was  correlated  to  the  slope  velocity  and  the  peak  velocity  as  meas¬ 
ured  on  the  two  highest  intensity  runs.  This  correlation  is  shown  in  Fig.  4.12. 
Using  this  correlation  peak  velocities  have  been  determined  from  the  distance  at 
the  fixed  time  for  all  laser  intensities.  These  peak  velocities  have  been  plot¬ 
ted  versus  peak  intensities  from  reference  1  in  Fig.  4.8  and  compared  with 
theory. 
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ABSTRACT 

A  simplified  model  is  developed  to  predict  the  impulse  imparted  to  a 
surface  when  irradiated  by  a  high  power  laser.  A  high  pressure  plasma  is 
created  by  a  laser  supported  detonation  (LSD)  wave  traveling  toward  the  energy 
source.  The  multi -dimensional  expansion  of  the  high  pressure  region  is 
described  by  patching  one,  two  and  three-dimensional  blast  wave  solutions  to 
the  LSD  wave.  The  patching  times  are  obtained  by  conserving  energy  between 
the  wave  solutions  and  are  shown  to  be  different  from  those  used  by  previous 
authors.  The  pressure-time-space  results  are  integrated  to  determine  the  net 
impulse  delivered  to  the  surface.  A  simplified  procedure  is  outlined  for  con- 
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verting  target  failure  criterion  to  laser  parameter  requirements 


1 


INTRODUCTION 


When  a  single  pulse  high  power  laser  irradiates  a  target  surface,  a 
laser  supported  detonation  wave  (LSD)  forms  above  the  target  surface  and 
propagates  into  the  background  gas.  The  high  pressure  gas  behind  the  LSD  wave 
transfers  momentum  to  the  target  and  represents  a  potential  damage  mechanism. 

A  measure  of  the  efficiency  with  which  momentum  is  imparted  to  the  target  is 
the  coupling  coefficient,  defined  as  the  impulse  per  unit  laser  energy.  Sev¬ 
eral  authors d-4)  have  sought  to  evaluate  and  optimize  the  coupling  coeffi¬ 
cient  with  respect  to  the  laser  beam  diameter,  pulse  duration  and  energy. 

The  impulse  imparted  to  the  target  during  irradiation  is  a  small  frac¬ 
tion  of  the  total  impulse  delivered.  Hence,  an  accurate  description  of  the 
momentum  transfer  must  include  the  decay  of  the  LSD  wave  into  a  hydrodynamic 
blast  wave.  The  transition  from  an  LSD  wave  to  a  blast  wave  has  been  describ¬ 
ed^1  f3,4)  by  patching  the  blast  wave  and  LSD  wave  solutions  together  at  an  in¬ 
termediate  time.  Hydrocode  simulations^)  have  confirmed  this  approach  but  a" 
unique  choice  for  the  patch  time  was  not  established.  It  is  shown  herein  that 
various  choices  for  the  patch  time  results  in  variations  in  the  energy  in  the 
blast  wave.  Since  the  energy  in  the  blast  wave  must  be  equal  to  the  laser 
energy  deposited  in  the  LSD  wave,  the  patch  time  may  be  established  by  spec¬ 
ifying  the  energy  in  the  asymptotic  blast  wave.  Using  this  technique  to  de¬ 
termine  the  patch  time,  it  is  shown  that  the  laser  energy  required  to  deliver 
a  fixed  impulse  to  a  target  may,  in  some  limits,  be  a  factor  of  five  greater 
than  previously  assessed. 


2.  PATCHED  SOLUTIONS 

The  high  pressure  plasma  generated  by  the  LSD  wave  imparts  momentum  to 
the  target  until  a  relaxation  wave  penetrates  the  plasma  from  the  surrounding 
low  pressure  ambient  environment.  The  relaxation  may  occur  from  the  axial 
(laser  beam)  direction  in  time  rz  or  laterally  in  time  T2d*  The  time  scales 
t2  and  T2d»  and  the  pressure  decay  laws  between  these  time  scales,  control  the 
total  momentum  imparted  to  the  target  by  the  laser.  The  LSD  wave  is  patched 
to  conventional  one,  two  and  three-dimensional  blast  wave  solutions  at  the 
times,  t2  and  T2d*  The  patching  sequence  varies  with  the  relative  magnitude 
of  these  time  scales. 

If  the  axial  relaxation  time  t2  is  less  than  T2d  (Case  I)  the  LSD  wave 
will  propagate  at  constant  velocity  and  pressure  until  time  r2  when  the  wave 
forms  a  one-dimensional  blast  wave  (1DSW).  The  1DBW  will  propagate  to  time 
T2d»  at  which  time  it  will  form  a  three-dimensional  blast  wave  ( 3DBW ) .  These 
patched  solutions  are  normalized  to  pressure  ps  which  represents  the  pressure 
of  the  LSD  wave  on  the  surface  of  the  target.  For  Case  I, 

t  <  T  :  p  =  p  ,  .  LSD  ( 1  ) 

z  s 

T  2/3 

Tz  <  t  <  T2D  :  P  =  ps  (^)  ,  1DBW  (2) 

Tz  2/3  T  6/5 

t2D  <  t  :  P  =  Ps  (— )  (— )  •  3DBW  (3) 
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-Similarly,  for  T2d  <  tz  (case  II),  the  LSD  wave  forms  a  two-dimensional  blast 


wave  (.2DBW)  at  time  T2d  and  a  3DBW  at  time  tz.  The  solutions  for  the  pressure 


in  case  II  are 


t  <  T. 


*  "  PS  ' 


1DLSD  (4) 


(  2D\ 

T2D  <  '  <  Tz  ‘  *  ■  Fs  (— >  ■ 


T  <  t 
Z 


p  -  ps  (^)  (J) 


3DBW  ( 6 ) 


The  critical  element  in  utilizing  the  patched  solution  is  that  of 


specifying  tz  and  T2d»  Pirri'1^  required  that  axial  relaxation  occurred  when 


the  laser  pulse  terminated,  i.e.,  Tz  equals  the  duration  of  the  laser  pulse. 


td.  Ferriter^2)  argued  that  axial  relaxation  could  not  occur  until  an  acous- . 


tic  wave  traveled  from  the  LSD  wave  to  the  target  surface,  notifying  the  tar¬ 


get  that  the  laser  pulse  was  off.  Neglecting  the  finite  velocity  of  the  gas 


behind  the  LSD  wave,  Ferriter  obtained  tz  =  2.83  Tp  whereas  tz  is  shown  to  be 
3.23  Tp  when  the  finite  velocity  of  the  gas  is  included. Similar  arguments 


have  been  utilized  to  evaluate  T2d*  Since  the  velocity  of  sound  behind  an  LSD 


wave  is  of  the  Order  of  the  velocity  of  the  wave,  Vw,  T2d  is  of  the  orcer  of 


Rs/vw  where  Rg  is  the  radius  of  the  irradiated  spot. 


The  above  choices  for  Tz  and  T2d  do  not  insure  that  the  asymptotic 


description  of  the  blast  wave  is  valid.  The  blast  wave  solution  for  the  pres¬ 


sure  may  be  expressed  as  p(E,t)  where  E  is  the  deposition  energy.  Various 
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choices  for  Tz  and  Tjd  in  Eqs.  (2),  (3),  (5),  and  (6)  will  inadvertently  alter 
the  energy  in  the  blast  wave.  This  suggests  that  there  must  be  appropriate 
values  of  Tz  and  t2D  that  constrain  the  energy  in  the  blast  wave  to  be  equal 
to  that  deposited  by  the  LSD  wave.  Below,  the  values  of  the  patch  times  are 
determined  via  this  criterion.  In  some  cases,  the  patch  times  and  blast  wave 
solutions  are  distinctly  different  from  those  discussed  above. 
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3.  ASYMPTOTIC  EVALUATION  OF  THE  PATCHING  TIMES 
The  variables  in  the  patched  pressure  solutions  [Eqs.  1-6]  are  deter¬ 
mined  directly  from  the  LSD  and  blast  wave  solutions.  The  LSD  wave  is  ex¬ 
pressed  in  terras  of  the  laser  intensity  IQf  gas  density  p,  and  the  ratio  of 
specific  heats  y.  The  velocity  Vw  and  pressure  py,  of  the  LSD  wave, (5)  and  the 
pressure  ps  on  the  target ^ )  are 

Vw  -  vv»,/3  ' 
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respectively,  where 


I2(Y2-1)]1/3  , 


a„  = 


and 


a.  = 


a1  /(Y+1)  , 


2y 

Y+1  \(y-i ) 


(— ) 
'  2y  ' 


The  time  scale  Tz  may  be  obtained  from  Eq.  (2)  and  the  known idbw 
solution  for  the  pressure  at  the  point  of  the  energy  deposition 


where  the  constant  has  been  evaluated  for  y  ■  1.2  and  is  the  energy  per 
unit  area  deposited  by  the  laser 


E1  *  *oTp  • 

The  above  expression  for  p(E^,  t)  at  the  target  surface  is  identical  to  that 
given  by  Eq.  (2)  if  and  only  if 


t2  -  dtp 


(Case  1)  (7) 


where 


a  «  0.6  , 

9  W 


for  y  *  1.2. 

Note  that  the  axial  relaxation  time  is  less  than  the  laser  pulse  dur¬ 
ation.  This  is  clearly  incorrect  when  viewed  locally  at  time  Tp.  However, 
the  dominant  momentum  transfer  occurs  during  the  asymptotic  portion  of  the 
blast  wave  and  the  local  details  are  sacrificed  in  order  to  accurately  predict 
the  asymptotic  behavior. 

when  the  lateral  relaxation  time  is  short  compared  .to  the  axial  relax¬ 
ation  time,  a  cylindrical  blast  wave  forms  above  the  target  surface  while  the 
LSD  wave  continues  to  move  toward  the  laser.  The  energy  deposited  per  unit 
length  perpendicular  to  the  target  surface  is 


E_  «  I  ttR  /V 
2  o  s  w 


and  the  asymptotic  blast  wave  solution  is6 


•for  T  *  1.2.  Comparing  this  expression  for  p(E2,t)  to  Eq.  (5),  T2d  becomes 


T2D  •  "V^w  '  (Case  XI)  (8) 

which  is  again  slightly  shorter  than  one  would  assess  from  a  local  analysis. 

The  axial  relaxation  time  Tz#  as  given  by  Eq.  (7),  is  valid  only  for 
Tz  <  T2d*  whereas  the  lateral  relaxation  time  T2d*  as  9i-ven  by  Eq.  (8),  is 
valid  only  for  T2d  <  tz.  In  both  cases,  the  expansion  fan  propagates  through 
the  hot,  high  density  gas  with  a  velocity  approximately  equal  to  Vw.  The 
length  scales  corresponding  to  Tz  and  T2d  are  vwTp  an<3  Rs»  respectively. 
Hence,  the  relaxation  times  are  approximately  Tp  and  Rs/Vw »  respectively. 

If  the  lateral  relaxation  occurs  first,  the  axial  relaxation  tz  must 


be  greater  than  given  by  Eq.  (7)  because  the  axial  expansion  is  now  occurring 
in  a  colder,  lower  density  gas.  Similarly,  if  the  axial  relaxation  occurs 


first,  the  lateral  relaxation  time  t2d  must  be  greater  than  given  by  Eq.  (8) 
because  now  the  lateral  expansion  is  occurring  in  a  colder,  lower  density  gas 
To  evaluate  these  limits,  both  Eqs.  (3)  and  (6)  are  required  to  be  identical 
to  the  asymptotic  3DBW  solution.  This  requirement  yields 

R3  V2 

X2D  =  — ~t]  (Case  I)  (9) 

T  VJ 
p  w 

and 

2 

t  =  cV  T  /R  .  (Case  II)  (10) 

z  w  p  s 

Equations  (7)— (10)  represent  a  complete  set  of  time  scales  for  the 
patched  pressure  solutions  [Eqs.  (1  )  —  (6)-)  and  clearly  define  the  range  of 
validity  of  Case  I  and  Case  II  in  terms  of  the  known  variables  Tp,  Rs  and  Vw. 


s 


f* 


ninnviiniriirvv  iiHUWiviniUMUliyuiUll 


Case  I  (Tz  .<  T2d):  Tp  <  Rs/Vw 

Case  II  (t2d  <  Tz):  Rs/Vw  <  Tp 

The  solution  for  the  total  momentum  imparted  to  a  target  surface  may 
be  generated  through  spatial  and  temporal  integrations  of  the  patched  pressure 
solutions.  To  accomplish  this,  patched  solutions  for  the  spatial  extent  of 
the  blast  wave  must  first  be  developed. 
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4.  PATCHED  SOLUTIONS  FOR  THE  SPATIAL  EXTENT  OF  THE  BLAST  WAVES 


The  expressions  for  the  lateral  extent  of  the  high  pressure  zone  are 
obtained  by  again  patching  blast  wave  solutions.  For  Tz  <  t2d  (Case  I),  the 
lateral  coordinate  R  is  equal  to  the  spot  radius  Rg  during  the  LSD  wave  and 
the  1DBW,  followed  by  the  3DBW  solution  for  t  >  t2d.  Ibis  is  expressed  as 


t  <  t2d  :  R  «  Rs 


LSD,  1DBW 


T2D  <  1 


t  2/5 

:  R  -  m  R  (—  ) 

1  S  ~2D 


where  m-|  is  a  constant  of  order  unity.  Equation  (12)  is  identical  to  the  3DBW 
solution  if  and  only  if 


mi  *  1.0  . 


Similarly,  for  Case  II  (t2e  <  Tz)>  the  lateral  extent  of  the  high 
pressure  zone  is  given  by 


t  <  T- 


R  =  Re  t 


LSD  (13) 


t  ' 

_  <  t  <  T  :  R  =  m_R  f - 1 


2DBW  (14) 


Matching  the  above  blast  wave  solutions  to  the  conventional  solutions®  yields 


m„  -  1.0 

2  • 


®3  “  1*0. 


The  2DBW  solution  appropriate  to  Case  II  is  valid  if  and  only  if  the 


two-dimensional  shock  radius  R  is  less  than  or  equal  to  the  target  radius  Rt. 


Beyond  Rt,  the  blast  wave  becomes  three-dimensional,  independent  of  TZJ  De¬ 


fining  Te  as  the  time  at  which  R  ■  Rt#  Eq.  (14)  yields 


T  «  aR  /R  V  , 
e  tr  s  w 


and  the  Case  II  solutions  are  valid  only  if  Tz  <  Te.  In  this  limit,  the  blast 


wave  undergoes  transition  from  2D  to  3D  before  it  reaches  the  edge  of  the 


target. 


Whenever  xe  <  tz,  the  blast  wave  undergoes  transition  from  2D  to  3D  at 


the  edge  of  the  target,  prior  to  time  tz.  in  this  limit,  modifications  to  the 


Case  II  solutions  are  necessary.  The  3DBW  solution  must  be  continuous  with 


Eqs.  (5)  and  (14)  at  R  *  R^.  The  3D  solutions  become 


T  6/5 

p  -  p.C^Hr) 
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3DBW  (17) 


T  1/2  2/5 
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2D  e 


3DBW  (18) 


respectively.  This  solution  does  not  allow  the  laser  and  blast  wave  energies 


to  be  identical.  Laser  energy  is  still  being  deposited  by  the  LSD  wave  when 


£1 

3 

S: 


k- 


K 


-the  cylindrical  blast  wave  reaches  the  edge  of  the  target.  Only  that  portion 
of  the  laser  energy  deposited  in  time  Te  is  contained  in  the  3DBW. 

The  above  limit  suggests  that  Case  II  can  be  restricted  to  T2d  <  xz 
<  Te  and  Case  III  be  defined  to  describe  the  limit  of  Te  <  Tz.  Expressions 
for  Tz  and  T2d  appropriate  to  Case  III  are  identical  to  those  in  Case  II.  The 
Case  III  solutions  for  p(t)  and  R(t)  are  the  Case  II  solutions/  with  Egs.  (17) 
and  (18)  replacing  Eqs.  (6)  and  (15).  respectively.  The  time  at  which  the 
2DBW  solution  transforms  into  a  3DBW  is  rfow  Te  rather  than  tz.  The  Case  I, 

II,  and  III  solutions  are  illustrated  in  Figs.  1-3,  respectively,  and  have 
been  expressed  in  a  form  that  will  be  used  to  determine  the  net  momentum 
transferred  to  the  target. 


5.  MOMENTUM  IMPARTED  TO  THE  TARGET 
The  total  impulse  (I)  imparted  to  the  target  is  the  spatial  and  tem¬ 
poral  integral  of  the  pressure  times  the  area. 

R 

I  *  /  /  p( t) 2nrdrdt  .  (19)  ' 

all  t  o 

The  patched  solutions  for  p(t)  and  R(t)  are  illustrated  in  Figs.  1-3  for  laaer 
pulses  of  increasing  time  duration.  The  Case  I  solutions,  appropriate  for 
short  laser  pulses,  are  illustrated  in  Fig.  1  and  depict  an  LSD  wave  trans¬ 
forming  into  a  1DBW  at  t2  and  then  into  a  3DBW  at  T2d*  For  intermediate 
pulses  (Fig.  2),  the  LSD  wave  transforms  into  a  2DBW  and  then  into  a  3DBW  be¬ 
fore  reaching  the  edge  of  the  target.  The  long  pulse  solution  (Fig.  3)  is 
identical  to  the  intermediate  pulse,  except  the  2DBW  reaches  the  edge  of  the 
target  before  decaying  into  a  3DBW.  Spatial  integration  of  the  pressure  over 
the  target  surface  must  be  terminated  at  either  the  laser  spot  radius  Rg, 
blast  wave  radius  R,  or  the  target  radius  Rt,  whichever  is  appropriate  at  time 
t.  Time  integrations  must  be  terminated  at  the  value  of  time  for  which  p(t) 
decreases  to  the  ambient  value,  p«,.  There  are  eight  distinct  subcases  corre¬ 
sponding  to  different  values  of  the  ambient  pressure  in  each  of  the  three 
cases.  They  are 

Case  la  :  P  +  P®  in  1DBW, 

Case  lb  :  p  p®  in  3DBW  with  R  <  Rt, 

Case  Ic  :  P  +  P«  in  3DBW  with*  R  >  Rt, 

Case  XXa  :  p  >  p.  in  2DBW, 

Case  lib  :  p  p®  in  3DBW  with  R  <  Rt, 

Case  I  Ic  s  p  p®  in  3DBW  with  R  >  Rt, 


Case  Ilia:  p  ♦  p»  in  2DBW, 

and 

Case  Illb:  p  *■  p.  in  3DBW  . 

Each  subcase  corresponds  to  different  values  of  the  laser  and  target 
parameters.  The  laser  parameters  specify  the  target  surface  pressure  ps,  wave 
velocity  Vw,  pulse  duration  Tp  and  spot  size  Rg.  The  regions  of  parameter 
space  corresponding  to  each  subcase  are  illustrated  in  Fig.  4.  Before 
elaborating  on  the  solutions  corresponding  to  each  of  the  eight  subcases,  it 
is  appropriate  to  discuss  the  approximations  made  in  the  spatial  and  temporal 
integrations  of  Eq.  (19). 

The  pressure  immediately  behind  the  blast  wave,  Pv,,  is  generally  a 
factor  of  three  greater  than  that  at  the  target  surface,  ps.  Since  the  pres¬ 
sure  decays  quite  rapidly  behind  the  blast  wave,  the  high  pressure  region  of 
the  pulse  is  not  included  in  the  spatial  integration  of  Eq.  (19).  This  ap¬ 
proximation  is  extremely  accurate  when  the  blast  wave  expands  beyond  the  tar-  ' 
get  and  only  the  uniform  pressure,  ps,  is  applied  to  the  target.  However,  the 
impulse  can  be  underestimated  by  "order  one"  if  the  blast  wave  does  not  expand 
beyond  the  target  (cases  lb,  Ila,  lib.  Ilia]. 

The  temporal  integration  of  Eq.  (19)  is  terminated  when  the  pressure 
on  the  target  surface  decays  to  the  ambient  value.  However,  the  blast  wave 
solution  is  valid  only  for  strong  shocks  (ps  >>  p»)  and  underestimates  the 
pressure  as  ps  ■+■  p«>.  This  results  in  an  underestimate  of  the  impulse  which  is 
appreciable  ("order  one")  only  if  the  integral  of  prdr  is  dominated  by  the 
limit  at  p».  This  occurs  only  when  the  blast  wave  does  not  expand  beyond  the 
target  surface  as  in  the  cases  cited  above. 


The  impulse  calculations  should  also  include  the  negating  effect  of 


the  ambient  pressure  opposing  that  of  ps  on  the  opposite  side  of  the  tar¬ 
get.  This  is  appreciable  (again,  "order  one")  when  the  integral  of  prdr  is 
dominated  by  the  limit  at  p»,  i.e.,  for  the  same  cases  as  cited  above.  While 
this  tends  to  cancel  the  two  effects  considered  above,  one  should  expect  as 
much  as  50%  errors  for  these  cases . 

Recognizing  the  approximations  involved,  the  eight  solutions  for  the 
total  impulse  imparted  to  the  target  are  given  below,  it  is  straightforward 
to  calculate  ps/p«o,  Rt/Rs»  vwTp/Rs  and  utilize  Fig.  4  to  determine  which  of 
the  following  solutions  is  appropriate. 

Defining 


p  tR  r  a  ' 
s  s  P 


(20) 


and 


R  3/2  R  1/2 

(tt-) 

s  w  p 


R.  3  R 
;  i.v  t  ; 

S  W  p 


1/6 


Vp- 


(21  ) 


the  solutions  are 


P  1/2 


Case  la: 


‘  ’  3 


-  2 


(22) 


Case  lb: 


c  P  1/2  .  R  1/2 

Srs’!  '  4  f  S'* 

c  -  )  +  -3(7*7")  '  2 


(23) 


w  p 


Case  Ic: 


a  \  l/2 

C*  +  -  2 


(24) 


P  R 
f  sy  s 


Case  Ila:  C  -  (— 

*•  w  p 


.  p  1/2  _  V  T 

Case  lib:  ,  -  -$(-**) 


Case  lie:  5  «  C*  -  j(  *-■**) 


Case  Ilia:  Same  as  Ila 


Case  Illb:  C  =  |  - 


w  p  s 


(VRs>' 
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P  /P 
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The  various  solutions  tabulated  above  illustrate  that  there  is  a  vari¬ 


ety  of  functional  relationships  between  impulse,  pulse  duration,  wave  veloc¬ 


ity,  spot  radius  and  target  radius.  Although  the  solution  in  each  regime  is 


relatively  simple,  the  boundaries  between  the  regimes  are  also  variable  and  it 


is  impossible  to  determine  a  general  scaling  law.  A  simple  computer  program 


has  been  written  which  utilizes  Fig.  4  to  determine  the  appropriate  case  and 


then  calculates  impulse  using  the  corresponding  analytic  solution.  Results 


are  illustrated  in  Fig.  5  and  are  compared  to  the  data  of  Kozlova  et  al.^J 


The  theoretical  predictions  agree  with  data  if  one  assumes  that  50%  of  the 


laser  energy  is  absorbed  by  the  gas.  The  particular  point  to  emphasize  here 


is  the  confirmation  of  the  value  of  0.6  for  the  constant  a.  Impulse 


scales  as  1/a  for  all  cases.  Local  theories (2,4)  predict  that  the  value  of 


is  of  the  order  of  three  and  the  corresponding  impulse  predictions  (via  this 


model)  are  a  factor  of  five  higher.  This  would  be  consistent  with  the  data  if 


and  only  if  the  energy  conversion  efficiency  were  reduced  to  10%, 
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The  reader  should  be  cautioned  against  concluding  that  the  predictions 
of  Ref.  4  are  always  a  factor  of  five  higher  than  those  obtained  here.  There 
are  several  differences  between  the  spatial  and  temporal  solutions  of  the  two 
models  and  there  is  only  one  limit  (case  la  with  Rg  *  Rt)  in  which  this  factor 
of  five  is  completely  apparent. 


6.  TARGET  DAMAGE  OPTIMIZATION 


The  analytic  solutions  for  the  impulse  imparted  by  a  high  energy  laser 
make  it  relatively  easy  to  study  the  effects  of  parametric  variation.  Con¬ 
sider  the  laser  energy  Eq,  pulse  duration  Tp  and  spot  radius  Rg  as  the  inde¬ 
pendent  laser  parameters.  The  beam  intensity  I0  is  expressed  as 
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and  the  wave  velocity  and  target  surface  pressure  become 
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respectively.  The  total  impulse,  can,  for  every  case,  be  expressed  as 

I  =  I(E0,  Tp,  Rt,  Rs/Rt)  .  (29) 

Ferriter^)  noted  that  maximum  impulse  is  delivered  when  the  laser 
beam  diameter- matches  the  size  of  the  target  (Rs  =  Rt) .  Expressing  the  total 
impulse  imparted  to  a  particular  target  R^  in  the  form  of  Eq.  (29),  it  can  be 
shown  that  all  expressions  for  I  (Eqs.  (22)— (28) }  approach  a  maximum  as 
Rs  *  Rt*  However,  the  regimes  in  which  these  solutions  are  valid  also  vary 
with  Rs/Rt  4/1(3  this  variation  is  not  considered  in  the  comparison.  The  re¬ 
sults  of  the  complete  model  are  illustrated  in  Fig.  6  for  a  104  joule  laser 


irradiating  a  10  cm  target.  These  calculations  clearly  demonstrate  that  maxi¬ 
mum  impulse  is  imparted  when  the  target  is  flood  loaded  (Rs  -  Rt) .  Such  a 
conclusion  is,  of  course,  subject  to  the  constraint  that  Rg  is  not  so  large  as 
to  reduce  the  beam  intensity  below  the  LSD  threshold,  nor  is  Rg  so  small  that 
the  plasma  is  overdense. 

In  the  limit  of  Rg  •  Rt,  Fig.  4  reduces  to  three  regions  corresponding 
to  Cases  la,  Ic  and  Illb.  The  expressions  for  the  impulse  per  unit  area  under 
the  spot  are  obtained  by  rewriting  Bqs.  (22),  (24),  and  (28),  retaining  only 
the  leading  term: 
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Structural  failure  of  a  target  requires  that  the  imparted  impulse  I/A 
exceeds  some  minimum  value  corresponding  to  some  minimum  pressure  ps.  A  typi¬ 
cal  failure  curve  is  illustrated  in  Fig.  7.  All  combinations  of  ps  and  I/A  in 
excess  of  the  critical  level  (curve  ABCD)  will  insure  failure.  Since  I/A  and 
ps  are  known  functions  of  (Ec,  Tp,  Rs),  it  is  relatively  straightforward  to 
invert  the  solutions  to  determine  the  E©  and  Tp  required  to  yield  a  specific 
combination  of  I/A  and  ps  on  spot  radius  Rg.  The  inversions  yield 
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Cases  la  and  Ic  correspond  to  regions  in  (I/A,  ps)  space  separated  by 


j  64oRs(pe#pa3) 


whereas  Case  Illb  reduces  to  a  single  curve  in  (I/A,  ps)  space  corresponding 
to  particular  values  of  Eo/Tp.  The  regions  are  illustrated  in  Fig.  7.  Param¬ 
eter  space  (ps,  I/A)  is  separated  into  the  three  regions  by  direct  application 
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of  Eqs.  (35)  and  (36).  That  portion  of  the  failure  curve  denoted  by  AB  cor¬ 
responds  to  Case  la.  This  portion  of  the  failure  curve  may  be  transformed 
into  (E0,  Tp)  space  via  Egs.  (30)  and  (31).  Similarly,  curve  BC  corresponds 
to  Case  Ic  and  may  be  transformed  into  (Ec,  Tp)  space  via  Eqs.  (32)  and  (33). 
The  single  point  on  the  failure  curve  denoted  by  C  corresponds  to  a  particular 
value  of  Eo/Tp  (Eq.  (34)]  and  completes  transformation  of  the  failure  curve 
into  (E0,  Tp)  space.  [Note  that  there  is  no  binary  combination  of  laser 
energy  and  pulse  duration  that  can  achieve  the  conditions  on  curve  CD] . 

A  relatively  simple  computational  procedure  [Eqs.  (30) -(36)]  has  been 
outlined  to  transform  the  failure  curve  (p8  vs.  I/A)  for  an  arbitrary  target 
into  the  corresponding  laser  requirements  (Eq  vs.  Tp) .  An  illustration  of 
this  procedure  is  presented  in  Fig.  8.  The  required  laser  energy  always  in¬ 
creases  with  increasing  pulse  duration.  This  is  a  consequence  of  the  fact 
that  the  imparted  impulse  always  decreases  with  increasing  pulse  duration. 

The  laser  failure  criterion  has  three  distinct  segments  corresponding  to  the 
three  possible  cases.  [Note  that  Case  la  need  not  exist  if  the  minimum  value 
of  I/A  for  target  failure  lies  in  the  Case  Ic  region) .  The  scaling  of  the 
laser  energy  with  the  constant  a  varies  in  each  region.  The  energy  scaling 
law  for  Case  Illb  is  a  function  of  the  shape  of  the  target  failure  curve.  The 
energy  scales  as  1/a  and  1/a^  in  Cases  la  and  Ic,  respectively.  The  asymptot¬ 
ically  derived  value  of  a(a  *  0.6)  is  a  factor  of  five  smaller  than  that  de¬ 
termined  from  local  analysis  (a  ■  3).  Hence,  the  laser  energy  required  to 
induce  structural  failure  of  the  target  may  be  five  or  more  times  greater  than 
that  which  would  be  determined  using  other  models.  Although  the  asymptotic 
value  of  a  «  0.6  was  confirmed  above  with  momentum  transfer  data,^)  its  im¬ 


plications  with  respect  to  the  required  laser  energy  have  not  been  confirmed. 
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7.  CONCLUSIONS 


An  existing  technique ( 1 “4 )  of  patching  blast  wave  solutions  together 
to  describe  the  pressure  field  above  a  laser  irradiated  target  has  been  modi¬ 
fied.  The  impulse  imparted  to  the  target  and  the  laser  energy  required  to  im¬ 
part  that  impulse  is  a  strong  function  of  the  patching  times .  Local  analyses 
do  not  uniquely  determine  these  times*  In  the  present  work,  the  conservation 
of  energy  is  used  to  determine  the  patching  times  in  a  more  rigorous  and  self- 
consistent  manner  (the  value  of  a  varies  20%  depending  upon  which  property  of 
the  blast  wave  is  matched).  Results  suggest  lower  impulse  per  unit  energy 
than  previously  suggested.  The  solutions  for  impulse  have  been  reduced  to 
simple  analytic  expressions  which  readily  invert  target  failure  criterion 
(pressure-impulse)  into  laser  requirements  (energy,  pulse  duration,  spot  size) 


for  target  failure 
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Fig.  1  Solutions  for  short  laser  pulses  (Tp  <  Rg /Vw,  Case  I). 


Fig.  2  Solutions  for  intermediate  laser  pulses  (Rs/Vw  <  Tp  <  Rt/Vw *  Case  II). 


Fig.  3  Solutions  for  long  laser  pulses  (Tp  >  Rt/Vw,  Case  III). 


Fig.  4  Range  of  parameter  space  appropriate  to  each  solution. 


Fig.  5  Comparison  of  predicted  impulse  to  data"?  for  a  «  0.6. 


Fig.  6  Impulse  imparted  by  a  104  joule  laser  on  10  cm  target. 


Fig.  7  Target  failure  criterion  (Rs  «  Rt) . 


Fig.  8  Laser  failure  criterion  (Rs  *  R^) . 
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